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I.  SUMMARY 


A  preliminary  analysis  of  the  vertical  motion  of  high  altitude 
balloons  was  presented  to  the  Office  of  Naval  Research  on  February  27, 
1961.  That  report  (Reference  1)  was  Technical  Report  I  of  a  series  of 
reports  which  have  been  sponsored  by  ONR  contract  Nonr-3164(00) . 

In  December  1963,  Technical  Report  II  (Reference  2)  was  published. 
This  report  extended  the  analysis  of  Reference  1  and  presented  a  computer 
program  which  describes  balloon  motion.  Technical  Report  III  which 
was  concerned  with  the  rotational  motion  of  high  altitude  balloons  was 
also  published  in  December  1963. 

Since  Technical  Report  II  was  published,  more  information  has 
become  available  for  correlation  purposes:  the  Stratoscope  series  has 
continued;  the  Fort  Churchill  summer  programs  have  produced  flight 
data  for  extreme  altitudes;  and  a  flight  was  made  by  ADL  with  the  help 
of  NCAR  which  measured  helium  gas  temperatures.  During  this  period  the 
analysis  and  the  computer  program  have  been  extended  and  improved.  They 
are  presented  in  this  report,  Technical  Report  IV. 

The  mathematical  formulation  of  the  model  is  presented  in  Section 
II. 


In  Section  III,  the  computer  program  is  described  in  great  detail. 
Since  this  version  of  the  program  is  more  accurate  than  that  presented 
in  Technical  Report  II,  all  program  users  are  hereby  requested  to  obtain 
this  newer  version.  With  the  permission  of  ONR,  the  program  is  available 
to  any  interested  party  in  the  subject  of  high  altitude  balloon  performance 

Finally,  in  Section  IV,  the  correlation  of  some  computed  and  actual 
flights  is  presented.  The  correlation  is  very  good.  The  results  from 
the  above  mentioned  ADL-NCAR  flight  have  increased  the  confidence  in 
the  validity  and  accuracy  of  the  computer  model. 


1 


3rtlmr  2).TUttle,Knr. 


II.  FORMULATION  OF  THE  THEORY  OF  BALLOON  DYNAMICS 


A.  INTRODUCTION 

Balloons  move  upwards  because  of  the  force  of  buoyancy,  a  simple 
and  well  understood  force.  The  upward  accelerating  force  (free  lift)  is 
equal  to  the  weight  of  the  displaced  air  minus  the  weight  of  gas  (gross 
list),  minus  the  total  weight  of  the  balloon  fabric,  payload,  ballast, 
etc.  So  it  is  a  simple  matter  to  write  the  dynamic  equation  of  vertical 
motion  of  a  balloon,  even  though  the  aerodynamic  effects  (induced  mass 
and  drag)  of  the  surrounding  air  introduce  some  complications. 

The  gross  lift,  however,  is  a  function  of  not  only  the  atmospheric 
conditions,  which  change  as  the  balloon  moves,  but  also  the  temperature 
of  the  gas,  which  depends  on  the  net  heat  received  by  the  balloon.  Due 
to  their  large  surface,  balloons  are  very  sensitive  to  heat  transfer, 
espet ialiy  to  thermal  radiation.  Thus,  a  balloon  acts  somewhat  like 
an  engine,  converting  heat  to  mechanical  energy,  and  its  equation  of 
motion  must  be  coupled  to  appropriate  energy  equations.  Writing  the 
appropriate  energy  equations,  which  will  account  for  all  important 
effects,  is  not  a  simple  matter. 

There  is  an  energy  equation  for  the  helium  gas.  The  gas  exchanges 
heat  with  the  fabric  through  free  convection  and  does  work  against  the 
atmosphere  when  it  expands.  The  expansion  of  the  gas  increases  the  lift 
and,  thus,  influences  the  motion  of  the  balloon.  This  equation  of  energy 
becomes  more  complicated  when  gas  is  valved  out  or  exhausted. 

There  is,  also,  an  equation  of  energy  for  the  fabric.  In  addition 
to  exchanging  heat  with  the  gas  through  free  convection,  the  fabric  ex¬ 
changes  heat  with  the  air  through  free  and  forced  convection.  Also,  the 
fabric  emits  thermal  radiation  and  receives  thermal  radiation  from  the 
sun,  the  earth,  and  the  atmosphere.  Because  of  the  importance  of  these 
effects,  the  various  radiation  fields  as  well  as  the  radiative  parameters 
of  the  fabric  must  be  specified  accurately. 

These  and  all  other  important  physical  processes  entering  balloon 
dynamics  are  described  in  detail  in  the  following  sections  and  are 
formulated  in  a  complete  system  of  equations.  A  few  words  about  the 
system  of  units  used  on  the  mathematical  formulation  of  the  problem  are 
in  order  here.  Length  is  in  feet  (ft),  time  in  seconds  (sec),  mass  in 
pounds  (lb),  temperature  in  degrees  Rankine  (  R) ,  and  force  in  pounds. 

All  other  quantities  are  expressed  in  terms  of  these  five  units.  For 
instance,  work  and  energy,  including  heat,  are  in  foot-pounds  (ft  lb). 

In  any  case,  the  units  of  every  quantity  are  given  in  the  List  of  Symbols. 

B.  THE  EQUATION  OF  VERTICAL  MOTION 

Let  wq,  wp,  and  wg  be  the  weights  in  lbs  of  the  payload,  balloon 
fabric  and  balloon  gas,  respectively.  The  total  mass  (in  slugs)  of 
the  balloon  system  as  it  moves  through  the  atmosphere  is: 
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— —  (wo  +  W„  +  VL 
g  G  F  B 


+  S  PA  VB) 


the  last  term  being  the  apparent  additional  mass  of  the  system  due  to 
the  surrounding  air.  VB  is  the  volume  of  the  balloon  in  ft3,  pA  is  the 
density  of  air  in  lbs/ft^,  and  g  =  32.2  ft/sec^.  Cg  is  a  constant  whose 
value  depends  on  the  shape  of  the  balloon.  For  a  spherical  balloon,  Cg 
is  equal  to  0.5. 

Let  z  denote  the  altitude  of  the  balloon  in  ft  and  y  its  velocity 
in  ft/sec.  Then: 


y  =  z  (1) 

and  the  equation  of  motion  of  the  balloon  system  is: 

<WG  +  WF  +  WB  +  S  P A  V  =  71  '  “3£“  °I>  PA  A,yly  (2) 

where  a  dot  denotes  differentiation  with  respect  to  time,  and  the  free 
life  FL  is  given  by: 


11  "  °A  Vg  '  UB  ‘  ”G  ‘  ”F  (3) 

The  second  term  of  the  right  hand  side  of  Equation  2  is  the  drag 
exerted  on  the  balloon  by  the  surrounding  air.  is  the  drag  coefficient 

and  A  is  the  effective  area  of  the  balloon  (in  ft ^)  in  the  direction  of 
motion.  For  a  spherical  balloon,  Cp  is  equal  to  about  0.45  for  Reynolds 
numbers  in  the  range  of  500  to  2  x  10-’  . 

In  balloon  terminology,  the  first  two  terms  of  the  right  hand  side 
of  Equation  3  are  known  as  gross  lift.  Using  the  equation  of  state,  the 
gross  lift  can  be  expressed  in  terms  of  quantities  which  are  well  known 
in  balloonry.  The  equation  of  state  (perfect  gas  law)  for  both  air  and 
gas  is: 

P  ■  PT  (A) 

2  3 

where  p  is  the  pressure  (in  lb/ft  ) ,  P  is  the  density  in  lb/ft  ,  T  is 
the  temperature  in  °R,  M  is  the  molecular  weight  and  R  is  the  universal 
gas  constant  in  ft  lb/lb  mol  R.  The  gross  lift,  GL,  can  now  be  expressed 
in  the  following  form: 
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GL 


(5) 


w„ 


,  ma  Pa<ta  +  e> 
*3  VPA+7T>J 


where  the  subscripts  A  and  B  pertain  to  air  ai yi  gas,  respectively.  The 
quantity  9  =  T^  -  T  is  known  as  the  gas  superheat  (actually,  it  is  a 
supertemperature) ,  while  ir  =  p  -  p  is  known  as  the  gas  superpressure. 
Since  the  percent  superheat  ana  superpressure  are  usually  small  quantities, 
Equation  5  can  be  reduced  to  the  following  approximate  relation: 


For  balloons  with  no  superpressure,  Equation  5  reduces  to: 


(5a) 


GL  =  w„ 


When  9  is  negative,  the  force  w 


Va 

is  negative  and,  in  balloon 


(5b) 


BMBTA 

terminology,  it  is  referred  to  as  thermodynamic  drag. 


C.  ENERGY  EQUATIONS 

.  For  a  closed  system,  tne  time  rate  of  change  of  its  internal  energy, 
U,  must  be  equal  to  the  rate  of  supply  of  energy  (in  heat  or  other  form), 
Q,  minus  the  rate  of  work  done  by  the  system,  W.  In  other  words: 

•  *  « 

U  =  Q  -  W  (6) 

The  units  in  this  equation  are  ft  lb/sec.  This  law  of  conservation  of 
energy  will  be  applied  to  the  balloon  gas  system  and  to  the  balloon 
fabric  system. 


For  the  balloon  gas  system; 


U  = 


dt 


C  w  T_ 
v  B  B 


C  w_  T  +  C  T_  w_ 
v  B  B  v  B  B 


(7) 


where  C  and  T_  are  the  specific  heat  at  constant  volume  in  ft  lb/lb  °R 
v  B  o 

and  the  temperature  in  R  of  the  gas,  respectively.  Notice  that  loss  of 
gas  by  exhausting  or  valving  is  included  in  the  above  equation. 

The  only  exchange  of  heat  of  the  gas  ic  with  the  balloon  fabric 
through  free  convection.  Let  q^  denote  the  rate  of  this  heat  transfer 
in  ft  lb/sec  from  the  gas  to  the  fabric.  In  addition  to  this  energy, 
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the  gas  system  also  loses  energy  at  the  rate  of  C  T  w  whenever  gas  is 
expelled  (w^  is  negative).  Therefore:  v  B  B 

Q  =  ~q6  +  Cv  TB  *B  <8> 

If  the  -ate  of  increase  of  the  balloon  gas  volume  is  V  ,  then  the 
rate  of  work  done  by  the  balloon  gas  on  the  atmosphere  is  pty  ,  where  p 
is  the  atmospheric  pressure  in  psf.  Also,  the  balloon  gas  does  work  on 
the  atmosphere,  whenever  gas  is  expelled,  at  the  rate  of  -pw^/pg,  where 
Pg  is  the  density  of  the  gas  in  lb/ft^.  Therefore: 


W  =  PVg  “  — p 


Now,  assuming  that  the  pressure  in  the  balloon  is  equal  to  the  atmospheric 
pressure  and  using  the  equation  of  state  (Equation  4)  for  the  balloon  gas. 
Equation  9  can  be  writtep  as: 


W  =  pvB  - 


T  w 
B  b 


Finally,  substituting  foi  U,  Q  and  W  from  Equations  7,  8  and  9a  in 
Equation  6,  the  following  energy  equation  for  the  gas  in  the  balloon  is 
obtained : 


Cv  "B  TB  *  -  I*  -  "VB  +  5.7  TB  “B 


For  the  balloon  fabric  system: 


4t-  (cf"ftf)  =  cf”fti 


where  C„  and  Tf  are  the  specific  heat  in  ft  lb/lb  R  and  the  temperature 
in  R  of  the  fabric,  respectively.  W  is  equal  to  zero  and  Q  is  given  by: 

Q  =  <12  "  q3  +  q4  f  q5  +  q6  +  q7  <12) 

where  all  the  q's  are  rates  of  heat  transfer  in  ft  lb/sec,  accounting 
for  the  effects  described  below. 


q^:  absorption  of  infrared  radiation  (from  earth)  by  fabric. 

q_:  emission  of  radiation  by  fabric. 
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q^:  forced  convection  to  fabric  from  air. 

q^:  free  convection  to  fabric  from  air. 

q^:  free  convection  to  fabric  from  balloon  gas. 

q^:  absorption  of  solar  radiation  by  fabric. 

Thus,  the  energy  equation  for  the  fabric  is: 


CF  WF  TF  'V  “  q3  +  q4  +  q5  +  q6  +  q7  C13) 

Equations  10  and  13  are  the  two  required  energy  equations.  Expressions 
for  the  q's  will  be  derived  in  Sections  G  and  H. 

D.  THE  EQUATION  OF  EXPANSION  OF  THE  BALLOON  GAS 

From  the  equation  of  state  (perfect  gas  law),  the  following 
expression  for  the  rate  of  expansion  of  the  balloon  gas,  V_,  is  obtained: 

D 

VB  '  ^  (”bTB+IbO  '  ~jf~  i  (14) 

It  should  be  kept  in  mind  that  the  pressure  of  the  balloon  gas  is  tal  m 
equal  to  the  atmospheric  pressure  p.  However,,  p  conforms  to  the  hydro¬ 
static  equation: 

•  • 

P  «-paZ 

Therefore,  substituting  for  p,  Equation  14  1 

VB  “  pMg  (WB  TB  +  TB  Wb)  + 

E.  VALVING,  EXHAUSTING  AND  BALLASTING 

The  pressure  in  inextensible  balloons  is  not  quite  the  same  as  the 
outside  atmospheric  pressure.  There  is  a  slight  overpressure  which 
causes  the  gas  to  flow  out  when  the  valve  is  open.  As  the  balloon 
reaches  its  ceiling  and  becomes  fully  inflated,  a  similar  small  pressure 
difference  across  the  appendix  of  the  balloon  causes  automatic  exhausting 
of  gas.  When  this  exhausting  is  sufficient,  it  can  prevent  the  balloon 
from  bursting  and  it  causes  stabilization  of  ceiling  without  high  alti¬ 
tude  bounce. 


lecomes 


A  v  • 
-  V_z 

p  B 


(15) 


(14a) 
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Let  E  and  V  be  the  rate  of  volumetric  gas  flow  in  ft3/sec  due  to 
exhausting  and  valving,  respectively.  Both  quantities  are  negative  for 
outflow  of  gas.  Then  the  rate  at  which  the  weight  of  the  balloon  gas 
changes  is  given  by: 

•  • 

*B  =  PB  (E  +  V)  <16> 

.  Valving  data  are  usually  given  as  lift  lost  per  unit  time  (lb/sec), 
Ly,  which  is  a  positive  quantity.  Mincer 

Lv  ”  -<PA  -PB>  V  (17) 


Equation  16  becomes: 


w„ 


(18) 


The  weight  of  the  payload,  wp,  is  changed  by  ballasting.  If  B  is 
the  ballasting  rate  in  lb/sec,  then: 


wG  -  -B  (19) 

F.  SPECIFICATION  OF  THE  ATMOSPHERE 

For  a  complete  specification  of  the  atmosphere  only  one  state 
variable  is  required  since  the  other  two  state  variables  can  be  computed 
from  the  hydrostatic  equation  and  the  equation  of  state.  Usually,  the 
atmosphere  is  specified  by  giving  a  temperature  profile  with  respect  to 
altitude.  The  pressure  and  density  can  then  be  computed  from  the  hydro¬ 
static  equation: 


= 

dz 


(20) 


and  the  equation  of  state: 


P 


P  A  TA 


(21) 


where  and  T*  are  the  molecular  weight  and  the  temperature  of  the 
atmosphere.  Eliminating  from  these  two  equations  and  integrating 
once,  one  obtains: 
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(22) 


A.  (\ 

R  J  ta 


z  A 
o 

where  p  Is  the  pressure  at  the  initial  altitude  z  .  It  is  clear  now 
that,  for  a  complete  specification  of  the  atmosphere,  the  initial  pressure 
p  (or  the  initial  density)  is  required  in  addition  to  the  temperature  profile 

Suppose  the  temperature  profile  is  specified  by  giving  a  number  of 
altitudes  and  corresponding  temperature  values,  z^  and  T  1  .  Then  in  the 
interval  n  and  n+1,  the  temperature  can  be  approximated  by  the  straight 
line: 


s  z  +  b 
n  n 


(23) 


where : 


T  , ,  -  T 
n+1  n 

z  ,  —  z 
n+1  n 


(24) 


Vn+l  -  Tn+lZn 

zn+l  '  2n 


(25) 


Integrating  Equation  22  with  T.  as  given  by  Equation  23,  the 
following  expressions  are  obtained.  For  an  isothermal  layer, 


and  s  ■  0: 
n 


£_ 

P„ 


exp 


MA 

■— T-  (*  "  Z) 

RT  n 
n 


(26) 


For  a  nonisothermal  layer,  Tn+^  J  and  s^  ^  0: 


£_ 

Pn 


(MA/snR) 


(27) 
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The  corresponding  expressions  for  the  density  are: 


P 


n 


exp 


(28) 


for  an  isothermal  layer,  and: 


(1  +  MA/snH> 


(29) 


for  a  npnisothermal  layer. 

G.  HEAT  TRANSFER  BY  CONVECTION 

1 .  Air-Side  Forced  Convection 


Forced  and  free  convection  as  well  as  radiation  depends  on  the 
shape  of  the  heated  object.  Convection  correlations  exist  in  the 
literature  for  plates,  cylinders  and  spheres  but  not  for  shapes  taken 
by  balloons,  which  can  be  anywhere  from  a  distorted  bubble  with  a  long 
stem  to  an  onion-like  shape  at  ceiling. 

The  balloon  is  considered  to  be  spherical  of  volume  Vg,  and  all 
heat  transfer  calculations  are  based  on  this  geometry.  Whenever  the 
deviation  of  the  actual  balloon  from  the  assumed  shape  is  thought  to 
have  an  important  effect  on  a  particular  heat  transfer  mechanism,  a 
correction  constant  is  introduced.  These  correction  constants  are 
evaluated  by  correlation  with  actual  flights. 

The  diameter,  D,  cross-sectional  area.  A,  and  surface  area,  S, 
of  a  sphere  of  volume  Vg  are  given  by: 


1/3 

D  =  1 . 24V 

u 

:•  (30) 

2/3 

A  -  1.21Vg 

(31) 

2/3 

S  =  4.83V_ 

D 

(32) 
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For  a  sphere  in  the  Reynolds  number  range  of  17  to  70,000, 
McAdams  (Ref.  3)  recommends  the  following  correlation  for  the  heat 
transfer  coefficient,  h,  by  forced  convection: 

—  =  0.37  (Re)  (33) 

where  k  is  the  conductivity  of  the  surrounding  medium  and  Re  is  the 
Reynolds  number.  He  points  out  that  turbulence  can  Increase  the  above 
value  of  h  by  40  to  60  percent.  The  laminar  flow  past  a  sphere  becomes 
turbulent  when  the  Reynolds  number  is  about  2.5  x  10^. 

For  all  big  balloons,  the  Reynolds  number  is  above  2.5  x  10** 
for  the  most  part  of  their  vertical  flight.  In  calculating  the  heat 
exchange  of  the  balloon  fabric  with  the  surrounding  air  by  forced  con¬ 
vection,  the.  balloon  is  assumed  to  be  spherical  and  Equation  33  is  used 
with  a  correction  constant  C^.  The  heat  transfer  is  given  by: 

q4  -  1.44  C4kAVB1/3  (Ta-Tf)(1.24  Vg173  y  ^  )0'6  (34) 

A 

where  kA  is  the  conductivity  of  air  in  ft  lb/ft  sec  °R  and  P  is  the 
viscosity  of  air  in  lb/.^t  sec.  A 

2.  .  Free  Convection 

For  free  convection  from  vertical  plates  and  inside  vertical 
cylinders,  McAdams  recommends  the  following  expression  for  the  heat 
transfer  coefficient,  h,  when  X  is  in  the  range  109  -  10. 


~  »  0.13  X173  (35) 

where  L  is  the  height  of  the  plate  or  cylinder.  X  is  the  product  of  the 
Grashof  and  Prandtl  numbers  and  it  is  proportional  to  the  third  power 
of  L. 

Let  us  assume  that,  with  an  appropriate  correction  constant, 
the  correlation  given  by  Equation  35  is  valid  inside  or  outside  a  sphere, 
as  long  as  X  based  on  the  diameter  of  the  sphere  is  within  or  near  the 
above  range. 

Taking  the  balloon  as  a  sphere,  it  can  be  shown  that,  for  both 
the  air-side  and  helium-side  of  large  balloons,  X  is  within  or  near  the 
above  range.  Using  Equation  35  with  correction  constants  C-  and  Cfi, 
the  heat  transfers  by  free  convection  in  the  air-side,  q^,  and  helium- 
side,  q^,  are  given  by: 
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1/3 


(36) 


0.628  Cr  k,  V 
5  A  B 


2/3 


(ta-tf) 
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0.628  C, 


kH  V] 


2/3 
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P  2  T 

H  B 


\/3 

PrH) 


(37) 


where  Pr^  and  Pr^  are  the  Prandtl  numbers  of  air  and  helium,  respectively, 
k^  is  the  conductivity  of  helium  in  ftlb/ftsec°R  and  is  its  viscosity 
in  lb/ftsec.  It  has  been  assumed  that  the  coefficients  of  expansion  of 
air  and  helium  are  equal  to  the  inverse  of  their  absolute  temperature. 

3.  The  Thermal  Parameters  of  Air  and  Helium 

In  the  following  paragraphs  a  brief  description  is  given  of  the 
Prandtl  number,  viscosity  and  thermal  conductivity  of  air  and  helium. 

The  temperature  range  of  interest  is  350  R  to  550  R. 

The  Prandtl  number  for  air  and  helium  is  essentially  constant 
and  equal  to  0.7  (Ref.  3). 

The  viscosity  is  a  function  of  the  temperature  only,  and,  in  a 
given  temperature  range,  it  can  be  taken  in  the  form  (Ref.  4) : 


v  =  AT 


(38) 


The  constants  A  and  n  are  determined  by  making  this  expression  conform  to 
two  experimentally  measured  values  of  the  viscosity  at  the  two  extreme 
temperatures  of  the  range  of  interest. 

-4 

From  Reference  5,  the  viscosity  of  air  is  1.333  x  10  and 
1.827  x  10"^  poise  at  temperatures  -69.4  and  18°C,  respectively.  Thus, 
Equation  38  becomes: 


PA  =  1.22  x  10‘6  Ta  0,883 


(39) 


where  the  viscosity  is  in  lb/ft  sec  and  the  temperature  in  R. 
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From  Reference  4,  the  viscosity  of  helium  is  equal  to  1.587x10 
and  1.967xl0_4  poise  at  temperatures  -60.9  and  17.6°C,  respectively. 
Therefore: 


H 


=  4.10  x  10 


0.682 


H 


(40) 


where  the  units  of  and  are  lb/ft  sec  and  °R,  respectively. 

The  specific  heat  at  constant  pressure  is  essentially  constant 
for  most  gases  (Ref.  4).  Therefore,  since  the  Prandtl  number  is  also 
constant,  the  thermal  conductivity  of  air  and  helium  must  depend  on 
temperature  in  the  same  way  as  the  viscosity,  i.e.,  in  the  form  of 
Equation  38  with  the  same  n.  The  only  unknown  constant  A  can  be 
determined  by  making  this  expression  conform  to  one  experimentally 
measured  value  of  the  thermal  conductivity  at  the  middle  of  the 
temperature  range  of  interest. 


From  Reference  4,  the  thermal  conductivity  of  air  ar.d  helium 


is  5.80  x  10" J  and 
ture  of  0°C.  Thus: 


3.52  x  10  4  cal/cm  sec  °K, 


respectively,  at  a  tempera- 


k.  =  4.08  x  10”7  T  0,883  (41) 

A  A 

ku  =  7.63  x  10'6  Tu  0-682  (42) 

H  rl 

where  the  thermal  conductivity  is  in  ft  lb/ft  sec  °R  and  the  temperature 
in  R. 


H.  HEAT  TRANSFER  BY  RADIATION 

1 .  Emission  of  Radiation  by  the  Balloon  Fabric 

The  fabric  emits  radiation  in  the  infrared  part  of  the  spectrum. 
The  flux  emitted  by  each  side  of  the  fabric  is  eoT^,  where  e  is  the 
average  emissivity  of  the  fabric  in  infrared  and  o  is  the  Stefan-Boltzmann 
constant.  Since  the  balloon  gas  is  considered  transparent,  part  of  the 
radiation  emitted  by  the  inner  side  of  the  fabric  is  absorbed  by  the 
fabric  so  that  the  effective  emissivity  e  -,  of  the  fabric  as  a  whole 
is  not  2e  but:  6 

2  3 

=  e  +  e  -  ea  -  eRa  -  eR  a  -  eR  a  -....  (43) 

where  a  and  R  are  the  average  absorptivity  and  reflectivity  of  the  fabric 
in  infrared,  respectively.  This  series  can  be  summed  to  give  the  following 
two  alternative  results: 
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(44) 


a 

T 

^  r 

__  £ 

1  1  r 

A  "  1  -  R 

r 

1  +  1  -  R 

r 

r 

where  the  subscript  r  is  used  to  denote  infrared  and  t  is  the  average 
transmissivity  of  the  balloon  fabric  in  the  infrared  spectrum  at  a 
reference  temperature. 


a  +  t  +  R  =1  (45) 

r  r  r 

If  the  balloon  is  considered  to  be  spherical,  the  radiating  fabric 
area  is  the  surface  of  a  sphere  of  volume  V_  as  given  by  Equation  32. 
Therefore,  the  effective  energy  emitted  by  the  fabric  is  given  by: 


q3  -  4.83  ero£  0VB  Tr4  (46) 

V 

where  o  is  in  ft  lb/ft  sec  °R  so  that  the  units  of  are  ft  lb/sec. 

Arthur  D.  Little,  Inc.,  has  developed  a  computer  program  for 
computing  the  effective  uultiple  absorptivity  and  emissivity  of  films. 
Using  measured  values  of  a  and  R  versus  wavelength  and  a  given  radiation 
intensity  spectrum,  1(1),  the  values  of  a  ,  e  and  R  of  the  film  are 
computed  from  expressions  such  as  the  folf owing:  r 

fad)  1(1)  dl 
ar  =  J*  1(1)  dl 


The  effective  absorptivity  and  emissivity  of  the  film  for  multiple  passes 
of  radiation  is  then  computed  from  Equation  44. 

This  program  is  used  to  calculate  the  effective  absorptivity 
and  emissivity  of  balloon  fabrics  in  infrared  as  well  as  their  effective 
absorptivity  in  the  sun’s  radiation  (0.3  to  3  microns).  For  infrared 
the  radiation  spectrum  used  is  that  of  a  black  body  at  a  temperature  of 
about  300°K.  For  solar  radiation, the  sun's  spectrum  outside  the  earth's 
atmosphere  is  used.  These  values  of  effective  absorptivity  and  emissivity 
are  used  for  the  entire  flight  of  the  balloon. 

There  are  two  sources  of  possible  errors  in  the  above  described 
method  of  computing  the  effective  absorptivity  and  emissivity  of  balloon 
fabrics.  Firstly,  the  spectra  of  the  various  radiations  are  not  the 
same  throughout  the  flight  of  a  balloon-  the  temperature  of  the  fabric 
changes,  and  the  spectra  of  the  infrared  radiation  from  the  earth  and 
of  the  solar  radiation  are  altered  through  the  atmosphere.  Secondly, 
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the  averaging  process  over  wavelength  should  be  carried  out  on  the 
multiple-pass  parameters  and  not  on  the  one-pass  parameters.  The  errors 
introduced  by  our  approximate  method  will  be  significant  when  the  radi¬ 
ative  properties  of  the  fabric  change  violently  with  wavelength  in  the 
spectrum  of  interest.  Indeed,  some  fabrics  behave  in  this  manner. 

2 .  Infrared  Radiation  From  the  Earth 

The  infrared  radiation  emitted  by  the  earth  is  absorbed  and 
emitted  again  by  the  atmosphere.  To  be  sure  the  resulting  radiation 
field  varies  throughout  the  atmosphere  not  only  in  its  total  intensity 
but  in  its  frequency  content  as  well.  Water  vapor  and  carbon  dioxide 
are  the  main  constituents  of  the  atmosphere  which  interact  heavily  with 
the  infrared  radiation  of  the  earth. 

There  are  nany  papers  dealing  with  the  absorbing  properties  of 
water  vapor  and  carbon  dioxide  layers.  For  instance,  Reference  6  gives 
extensive  tables  and  curves  describing  these  properties.  From  these 
data  it  is  possible  to  construct  a  theoretical  model  for. the  computation 
of  the  radiation  received  by  a  balloon  as  it  moves  through  the  atmosphere. 
However,  such  a  model  would  require  an  inhibitiv?  amount  of  computation 
and  it  would  be  almost  useless  since  the  concentration  and  stratification 
of  water  vapor  and  carbon  dioxide  in  the  atmosphere  is  quite  unpredictable. 

The  infrared  radiation  field  in  the  atmosphere  has  been  measured 
by  many  observers.  Most  important  in  this  field  are  the  works  of  Gergen 
(Refs.  7-11).  The  radiation  field  versus  altitude  is  measured  with  black 
ball  flights.  The  equilibrium  temperature  of  the  black  ball,  T  ,  is 
recorded  and  given  in  charts  versus  altitude. 

This  type  of  measurement  has  been  made  for  many  locations.  The 
results  show  that  the  profile  of  T  varies  with  geographical  and  seasonal 
conditions.  Essential  changes  of  Ehe  T  profile  for  the  same  location 
can  take  place  within  days.  The  measurement  is  carried  out  during  the 
night,  so  that  the  detector  is  not  effected  by  the  sun.  It  has  been 
estimated  that  the  day  radiation  field  is  about  10  F  higher  than  the 
measured  night  field. 

The  results  of  these  measurements  support  the  following  approxi¬ 
mate  but  simple  general  rule.  At  ground,  T  is  usually  less  than  the 
temperature  of  the  air,  the  deviation  being  not  more  than  about  10  F. 

Thgn  T  dectaases  almost  linearly  with  altitude  up  to  tropopause,  where 
ctT  israbout  307,  of  its  value  at  ground.  From  there  to  higher  altitudes, 

T  remains  approximately  constant. 

3 .  Absorption  of  Infrared  Radiation  by  the  Balloon  Fabric 

The  use  of  a  spherical  detector  to  measure  the  infrared  radiation 
field  is,  indeed,  a  very  fortunate  coincidence.  Since  in  the  theoretical 
model  the  balloon  is  considered  to  be  spherical,  i.e.,  of  the  same  shape 
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as^the  detector,  the  radiation  flux  incident  on  the  balloon  is  simply 
0 ^ j.  and  the  energy  absorbed  by  the  fabric  can  be  calculated  very  easily. 
This  would  not  be  the  case, if  the  detector  had  a  different  shape.  In 
the  following  calculation  of  the  infrared  radiation  absorbed  by  the 
fabric,  it  is  assumed  that  the  T  profile  is  known.  If  this  profile 
has  not  been  measured  for  a  particular  location  and  time,  then  the  simple 
rule,  described  in  the  last  paragraph  of  the  preceding  section,  is 
recommended  for  its  specification. 

Since  the  balloon  gas  is  considered  transparent,  the  absorption 
of  the  incident  flux  oT^  by  the  fabric  is  multiple,  as  in  the  case  of 
emission  described  in  Section  1.  Thus,  an  effective  absorptivity, 

£,  and  not  the  one-way  absorptivity  of  the  fabric  must  be  used.  The 
cafculation  of  a  is  done  as  described  in  Section  1  with  a  radiation 
spectrum  identical  to  that  of  a  black  body  at  temperature  T  .  With  the 
absorbing  surface  S  as  given  by  Equation  32,  the  infrared  radiation 
energy  absorbed  by  the  fabric  is  given  by: 


q2 


4.83 


a  r  a  V„ 
ref  B 


2/3 


4.  Solar  Radiation 


(47) 


The  total  intensity  of  solar  radiation  outside  the  earth's 
atmosphere  is  2.05  cal/cm  min  (Ref.  6)  or  96  ft  lb/ft  sec.  The  peak  of 
i;-s  spectrum  is  in  the  visible  range  of  frequencies.  A  considerable 
amount  of  energy  is  in  the  infrared  and  some  energy  is  in  the  ultraviolet. 

As  this  energy  enters  the  atmosphere,  part  of  it  is  absorbed  by 
the  various  atmospheric  constituents.  Ozone  absorbs  ultraviolet;  ozone, 
oxygen  and  water  vapor  absorb  visible;  water  vapor  and  carbon  dioxide 
absorb  infrared.  By  the  time  the  solar  radiation  reaches  the  earth  its 
total  intensity  is  reduced  to  1.44  cal/cm^min,  when  the  azymuth  angle 
of  the  sun  is  zero  (i.e.,  the  sun  is  directly  overhead),  and  its  frequency 
spectrum  is  altered. 

An  approximate  way  of  describing  the  attenuation  of  the  solar 
radiation  in  the  atmosphere  is  through  the  optical  air  mass,  m,  which 
is  the  ratio  of  the  path  length  of  the  sun's  rays  through  the  atmosphere 
to  the  normal  path  length.  The  following  table  is  constructed  from  data 
given  in  Reference  6. 
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Table  1. 


Variation  of  Solar  Radiation  (Relative  to  Radiation 
Outside  Atmosphere)  With  Optical  Air  Mass. 


Optical  Air  Mass  (m) 

0 

1 

2 

3 

4 


Solar  Radiation  (T) 

1.0 

.7018 

.5596 

.4595 

.3849 

.3249 


For  a  computer  program,  it  is  expedient  to  have  closed  form 
expressions  instead  of  tables  of  values.  The  following  expression 
conforms  to  the  values  of  the  above  table  with  an  error  of  less  than  3%. 


T  -  1/2  (e“°'65m  +  e“°*095m^  (48) 

The  optical  air  mass,  m,  depends  on  the  altitude  and  on  the 
sun's  azymuth  angle,  6.  When  the  sun  is  directly  overhead,  its  azymuth 
angle  is  zero.  Sea  level  sunrise  or  sunset  correspond  to  6  =  90°.  For 
6  =  0,  m  is  equal  to  0  and  1  for  a  point  outside  the  atmosphere  and  at 
sea  level,  respectively. 

Let  mQ  be  the  value  of  m  at  sea  level.  Then,  keeping  6  constant, 
the  value  of  m  at  any  altitude,  where  the  pressure  is  p,  is  given  by: 


m  =  m 


°  P„ 


where  is  the  pressure  at  sea  level. 


(49) 


The  variation  of  m  with  5  is  given  in  Table  16-18  of 
Reference  6.  It  can  be  ver?fied  easily  that  the  following  expression 
conforms  to  the  values  of  this  table  with  an  error  of  less  than  0.5%. 


m 


1228.6  +  (613.8  cos  6  )‘ 


1/2 


-613.8  cos  o$6  £  90  (50) 


Thus,  for  given  time  (and,  therefore,  given 6  )  and  given 
altitude,  the  optical  air  mass  can  be  found  from  Equations  49  and  50. 
Then  the  solar  radiation  (relative  to  that  outside  the  atmosphere)  for 
that  particular  time  and  point  in  the  atmosphere  can  be  computed  from 
Equation  48.  Obviously,  this  procedure  is  valid  for  values  of  6  up  to 
90  .  Sea  level  sunrise  and  sunset  will  occur  at  6  =  90°,  and  a  point 
at  sea  level  will  not  receive  any  solar  radiation  for  6  greater  than 
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90  .  This  is  not  the  case  for  a  point  higher 
where  sunrise  and  suusat  occur  at  values  of  <$ 
point  will  receive  solar  radiation  for  values 
A  procedure  must  be  devised  for  computing  the 
range  of  6. 


up  in  the  atmosphere, 
greater  than  90°.  This 
of  6  greater  than  90°. 
solar  radiation  in  this 


Consider  point  A  in  the  atmosphere  (Figure  1),  which  has  an 
altitude  z.  At  this  particular  time,  the  sun's  rays  are  along  line  SA, 
and  6  is  larger  than  90  .  The  atmospheric  path  travelled  by  the  sun's 
rays  is  CA.  Line  BD,  drawn  from  the  earth's  center,  is  perpendicular 
to  line  AC.  The  altitude  of  point  B,  z^,  is  given  by: 


z^  =  (R  +  z)  sin  6  -  R  (51) 

where  R  is  the  radius  of  the  earth. 

For  point  B,  the  azymuth  angle  of  the  sun  is  equal  to  90°.  Thus, 
the  optical  air  mass  for  point  B,  m1  ,  can  be  computed  from  Equations  49 
and  50  with  -6  =  90°  and  p  =  p^,  where  p^  is  the  pressure  at  z  .  Then 
the  attenuation  of  solar  radiation  at  point  B  is  T(rn^)  as  give^by 
Equation  48. 


In  travelling  the  remaining  air  path  BA,  the  solar  radiation  is 
reduced  further.  Let  m„  be  the  optical  air  mass  corresponding  to  altitude 
z  and  azymuth  angle  of  the  sun  equal  to  180  -  6,  i.e.,  when  the  sun  is 
at  S'.  Any  radiation  originating  at  B  is  reduced  to  T(mL)  when  it  arrives 
at  C' ,  and  any  radiation  originating  at  S'  (or  C')  is  reduced  to  T^p 
when  it  arrives  at  A.  Therefore,  any  radiation  originating  at  B  is 
reduced  to  T(m  )/T(m  )  when  it  arrives  at  A.  Thus,  the  solar  radiation 
received  by  point  A  (relative  to  that  outside  the  atmosphere)  is  given 
by: 


T2(mL) 

T(«n2) 


(52) 


Tne  range  of  6  covered  by  this  computation  is: 


90°  <  6  90°  +  cos"1  (53) 

K  “r  Z 

For  larger  values  of  6  ,  the  solar  radiation  received  by  point  A  is 
equal  to  zero. 

5 .  Absorption  of  Solar  Radiation  By  the  Balloon  Fabric 

The  area  of  a  spherical  balloon  normal  to  the  solar  ray c,  which 
come  from  infinity,  is  equal  to  the  area  of  a  circle,  A,  as  given  by 
Equation  31.  This  is  the  effective  absorbing  area  of  the  fabric  for  one 
pass  of  radiation.  However,  since  the  balloon  gas  is  again  considered 
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transparent,  there  is  multiple  absorption  by  the  fabric  and  an  effective 
absorptivity  a  must  be  used,  not  the  one-way  absorptivity  of  the 
fabric.  Ihe  subscript  v  stands  for  solar  radiation.  This  absorptivity 
is  computed  as  described  in  Section  1  with  the  solar  radiation  spectrum. 

Thus,  the  solar  energy  absorbed  by  the  fabric  is  given  by: 


q  =116  a  Tv  2/3  (54) 

m7  vef  B 

The  numerical  constant  in  this  equation  comes  from  the  solar  constant 
outside  the  atmosphere  (96  ft  lb/ft2sec)  times  the  numerical  constant 
of  the  effective  area  (1.21)  as  given  by  Equation  31.  T  is  a  function 
of  the  optical  air  mass  computed  as  described  in  Section  4. 

It  remains  to  define  the  sun's  azymuth  angle,  6  ,  in  a  time 
coordinate  system  which  is  related  to  the  flight  of  the  balloon.  Let 
LONG  and  LAT  be  the  longitude  and  latitude  of  the  balloon  in  degrees, 
respectively,  If  at  time  zero,  the  Greenwich  hour  angle  is  GHA,  then 
at  any  subsequent  time  t  the  local  hour  angle,  LHA,  is  given  by. 

LHA  =  GHA  -  LONG  +  “  (55) 


where  LHA  and  GHA  are  in  degrees  and  t  is  in  seconds.  Let  DEC  be  the 
declination  sf  the  sun  in  degrees.  DEC  is  a  slowly  varying  function  of 
time  so  that  it  can  be  considered  constant  throughout  a  balloon  flight. 
Then  6  is  given  by  the  following  equation. 

cos 6  =  sin  LAT  sin  DEC  +  cos  LAT  cos  DEC  cos  LHA  (56) 


I.  SUMMARY  OF  EQUATIONS.  REQUIREMENTS  FOR  UNIQUE  SOLUTION^ 

In  summary,  the  governing  differential  equations  of  the  model  are: 

,  .  „  (57) 


(w„  +  w,  +  wB  +  CB  PA  VB)  *  =  PA  VB  -  wG  -  wF  -  wB 


~G  F 

1.21  „  II  v  2/3 
2g  D  Aiyiy  B 


Cv  WB  TB  =  TL  *B  -  PVB  '  q6 


(58) 

(59) 


CF  WF  TF  =  q2  ‘  q3  +  q4  +  q5  +  q6  +  q7 


(60) 
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VB  =  <UB  +  TB  *8)  +  T4  VB  * 


(61) 


w  =  p  R  -  -  L 

B  B  PA  '  PB  V 


(62) 


Notice  that  in  the  drag  term  of  Equation  58,  the  effective  area  has  been 
taken  equal  to  the  cross-sectional  area  of  a  spherical  balloon  of  volume 


Thus,  there  are  six  first  order  nonlinear  differential  equations. 

The  independent  variable  is,  of  course,  tim».  Let  the  six  dependent 
variables  be  z,  y,  w^,  Vg,  Tg,  and  Tp.  For  a  unique  solution  of  these 
variables,  all  other  variables  must  be  known.  In  detail,  the  following 
items  must  be  specified. 

1.  The  launch  site  (longitude  and  latitude)  and  the  launch  time. 

2.  The  atmospheric  temperature  (T^) ,  pressure  (p),  density  (P  ), 
thermal  conductivity  (k^)  and  viscosity  ( M  ).  The  specifica¬ 
tion  of  all  these  properties  is  accomplished  by  means  of  a 
temperature  profile  (T^  vs  z)  and  the  initial  pressure  or 
density. 

3.  The  infrared  radiation  field  in  the  atmosphere.  It  is  specified 
by  giving  the  "black  ball  radiation  equilibrium  temperature" 
profile  (Tf  vs  z) . 

4.  The  weight  (wp  ,  specific  heat  (C^.) ,  infrared  effective  multiple 

emissivity  ( e  r  p  and  effective  multiple  absorptivity  for 
infrared  (a  rJ  and  for  solar  radiation  ( a  ,.)  of  the  balloon 
fabric.  ret  ve£ 

5.  The  total  payload  including  ballast  (w^,) ,  and  the  ballasting 
schedule. 

6.  The  valving  schedule  and  exhausting  rate.  Specification  of  the 
latter  is  very  difficult,  if  not  impossible,  since,  in  addition 
to  the  size  and  shape  of  the  appendix,  the  instantaneous  pressure 
differential  across  the  appendix  must  be  known.  In  the  computer 
program,  exhausting  is  allowed  through  a  mathematical  expedient 
which  will  be  described  in  Section  III-C. 

7.  The  thermal  conductivity  (k^)  and  viscosity  ( u ^)  of  the  gas 

as  functions  of  the  gas  temperature  (T,,)  . 

Jd 
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8.  The  Prandtl  number  for  air  and  for  the  balloon  gas.  Also,  the 
specific  heat  at  constant  volume  of  the  gas  (C  )  and  R,  M  and 

V  V  A 

9.  The  five  constants  Cg  C^,  C^,  C,.  and  Cg.  As  it  will  be  shown 
in  Section  IV  good  cbrrelation  with  actual  flights  is  obtained 
with  the  following  values  for  these  constants: 

S  ■  °-5 


C.  =  1.5 

4 


c6  -  1.0 

Furthermore,  for  a  unique  solution,  initial  values  of  the  six 
dependent  variables  must  be  given.  For  a  balloon  flight,  the  initial 
values  of  z,  y,  T  and  T_,  can  be  obtained  rather  easily  with  the 
exception,  perhaps,  of  T„.  As  for  the  initial  values  w  and  Vg,  they 
can  be  calculated  from  the  initial  value  oi  the  free  lift,  which  must 
be  given,  and  the  equation  of  state  with,  of  course,  the  gas  pressure 
equal  to  the  atmospheric  pressure. 
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III.  THE  COMPUTER  PROGRAM 


A.  INTRODUCTION 


A  computer  program  has  been  devised  to  solve  the  problem  of  the 
vertical  flight  of  a  balloon  as  formulated  in  Section  II.  The  program 
integrates  numerically  the  six  first  order  differential  equations 
(Equations  57  to  62)  once  the  necessary  input  data  are  given. 

The  program  is  coded  in  the  FORTRAN  II  algebraic  language  and  it 
has  been  operated  in  the  IBM  7090  Data  Processing  Equipment  in  con¬ 
junction  with  the  FORTRAN  MONITOR  SYSTEM.  The  computer  time  depends 
primarily  on  the  integration  time  interval,  the  on  and  off  line  print¬ 
ing  period  and  the  plotting  period.  For  an  integration  time  interval 
of  20  seconds,  an  off-line  printing  period  of  15  minutes  and  a  plotting 
period  of  5  minutes,  the  times  required  by  the  7090,  which  does  the  in¬ 
tegration,  and  by  the  1401,  which  prints  the  output,  are  both  about  one 
hundreth  Of  the  actual  flight  time. 

In  the  following  sections,  the  program  is  explained  in  detail. 

The  FORTRAN  listing  of  the  entire  program  and  a  general  flow  chart  are 
given  at  the  end  of  this  report.  The  symbols  used  in  the  program  are 
defined  in  the  List  of  Symbols.  The  flow  chart  is  intended  to  show  in 
one  compact  picture  the  general  f} ow  of  information,  sequence  of  opera¬ 
tions  and  logic  of  the  program. 

B .  INPUT  DATA 


The  input  data  contain  all  the  information  required  for  a  unique 
solution  (see  Section  II-I)  as  well  as  information  pertaining  to  the 
integration,  and  to  printing  and  plotting  of  output. 

All  input  data  are  read  in  the  program  from  the  input  tape,  which 
gets  this  information  from  a  deck  of  cards.  No  information  can  be  fed 
to  the  computer  on-line.  However,  all  six  sense  switches  can  be  used 
during  the  course  of  computation  to  perform  various  functions. 

Besides  the  leading  title  card,  which  contains  alphanumeric  infor¬ 
mation  in  format  12A6,  there  are  two  categories  of  input  cards:  control 
cards  and  data  cards.  Control  cards  must  have  only  an  integer  number, 
from  0  to  13,  in  the  first  two  columns  in  format  12.  Data  cards  must 
have  the  first  two  columns  blank  and  then,  in  successive  fields  of  10, 
input  data  in  format  F10.4.  The  data  are  divided  in  eleven  types. 

Each  type  must  be  preceded  by  a  control  card. 

The  content  of  a  complete  set  of  data  cards  for  one  flight  is 
shown  in  Table  2.  Each  line  represents  one  card.  Notice  the  order  of 
control  cards.  The  control  card  with  the  number  1  will  make  the 
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computer  to  start  computing  and,  therefore,  it  must  be  placed  after  all 
the  data  cards.  The  control  card  with  the  number  2  will  make  the  com¬ 
puter  to  CALL  EXIT  instantly  and,  therefore,  it  must  be  the  last  card 
in  the  deck.  It  is  not  necessary  to  have  the  various  types  of  data  in 
the  indicated  order.  They  can  be  rearranged  as  long  as  they  carry  the 
correct  control  card.  For  instance,  data  type  4  can  be  placed  after 
data  type  8,  and  so  on.  However,  within  a  given  type  of  data,  the  data 
must  be  given  in  the  indicated  order.  A  control  card  must  always  be 
followed  by  the  appropriate  data  cards. 

Type  3  data  contain  the  specific  heat  of  the  fabric  (Cl)  in 
ftlb/lb6 R,  the  constants  of  convective  heat  transfer  (C4,  C5  and  C6) 
which  are  equivalent  to  the  constants  C^,  C5  and  Ct  of  Section  II,  the 
infrared  effective  multiple  absorptivity  (ABIR)  and  emissivity  (EMIR) 
of  the  fabric,  and  the  effective  multiple  absorptivity  for  solar  radi¬ 
ation  (ABUV)  of  the  fabric. 

Type  4  data  contain  the  initial  weight  of  the  payload  including 
ballast  (WGO)  in  lbs,  the  weight  of  fabric  (WFO)  in  lbs,  die  inflated 
volume  of  the  balloon  (VBM)  in  ft3  as  specified  by  its  manufacturer, 
the  initial  free  lift  (FLO)  in  lbs,  and  the  aerodynamic  constants  for 
drag  (CD)  and  apparent  additional  mass  (CB)  which  are  equivalent  to  the 
constants  C^  and  of  Section  II. 

Type  5  data  contain  the  initial  temperature  ofQthe  gas  (TBO)  in 
°R,  the  initial  temperature  of  the  fabric  (TFO)  in  R,  the  initial 
altitude  of  the  balloon  above  sea  level  (ZO)  in  ft,  the  initial  velocity 
of  the  balloon  (Y20)  in  ft/sec,  the  declination  of  the  sun  (DEC)  in 
degrees,  and  the  latitude  (XL AT)  and  longitude  (XL0NG)  of  the  balloon 
in  degrees. 

Type  6  data  contain  the  time  interval  for  the  integration  (H)  in 
sec,  the  printing  period  of  output  (DPR)  in  sec,  the  initial  time  (XO) 
in  sec,  the  final  time  (XT)  in  sec,  the  Greenwich  mean  time  (GMTS)  in 
sec  at  the  initial  time,  the  „*eenwich  hour  angle  (GHA)  in  degrees  at 
the  initial  time,  and  the  plotting  period  of  output  (DP)  in  sec. 

Type  7  data  contain  two  arrays  of  time  (TCSI)  in  sec  and  correspond¬ 
ing  values  of  the  so  called  "solar  radiation  factor"  (CCSI) .  They  must 
be  in  order  of  increasing  TCSI.  The  number  of  cards  must  not  exceed 
100.  Solar  radiation  for  a  clear  sky  is  computed  automatically.  Re¬ 
duction  of  the  solar  radiation  caused  by  transient  effects  such  as 
clouds  are  accounted  for  by  the  solar  radiation  factor.  When  no  such 
effects  are  considered,  th <  correct  and  sufficient  data  of  this  type 
are  the  following  two  cards: 

TCSI(l)  =  XO  CCSI(l)  =1.0 

TCSI(2)  =  XT  CCSI(2)  =1.0 
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Table  2.  Content  and  Order  of  Input  Cards  for  One  Flight 


Field 

Format 

_ 

Ml... 

22| 

32|  42| 

52| 

62| 

72| 

TCSI(l)  CCSI(l) 


TCSI(NSl)  CCS I (NS I) 


VO 


TV(1)  W(l) 


TV(NVI)  W(NVI) 


BO 


TB(1)  BB(1) 


TB(NBI)  BB(NBI) 


EMIR  ABUV 


WGO 

WFO 

VBM 

FLO 

CD 

CB 

XL AT  XL0NG 


12A6 


12 


F10.4 


12 


F10.4 


12 


F10.4 


12 


F10,4 


12 


F10.4 


F10.4 


F10.4 


12 


F10.4 


F10.4 


F10.4 


F10.4 


12 


F10.4 


F10.4 


F10.4 


F10.4 


Field 


Format 
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To  account  for  a  transient  effect,  more  cards  are  needed  containing 
values  of  CCSI  smaller  than  1.0  and  the  corresponding  times  TCSI. 

Type  8  data  contain  the  "manual"*  valving  rate  (VO)  in  pounds  of 
lost  lift  per  second  on  the  first  card.  Even  when  no  gas  will  be  valved 
manually  in  the  course  of  computation,  this  card  must  lead  the  valving 
data.  The  following  cards  contain  two  arrays  of  time  (TV)  in  sec  and 
corresponding  values  of  automatic  valving  rate  (W)  in  pounds  of  lost 
lift  per  second.  They  must  be  in  order  of  increasing  TV.  The  number 
of  cards  must  not  exceed  100.  Also,  TV(1)  must  be  equal  to  or  smaller 
than  XO  and  TV(NVI)  must  be  equal  to  or  greater  than  XT.  The  following 
example  shows  the  meaning  of  these  valving  data  and  the  way  they  are 
used  by  the  program. 

Suppose  for  a  flight  XO  =  0  and  XT  =  36,000  sec,  and  that  the 
valving  data  for  this  flight  are: 


VO 

=  1.0 

TV(1) 

=  0 

W(l)  =  0. 

TV(2) 

=  18000. 

VV(2)  =0.3 

TV  (3) 

=  18100. 

3 

/-N 

U> 

II 

O 

• 

TV  (4) 

=  36000. 

VV(4)  =  0. 

These  data  mean  that  the  "manual"  valving  rate  is  1  and  that  automatic 
valving  takes  place  only  between  times  18000  and  18100  at  a  constant 
rate  of  0.3.  Thus,  the  total  amount  of  gas  valved  automatically  is 
equivalent  to  a  loss  of  30  lbs  of  lift. 

Type  9  data  contain  ballasting  information.  The  first  card  con¬ 
tains  the  "manual"*  ballasting  rate  (BO)  in  lbs/sec.  Even  when  no 
manual  ballasting  will  be  done  in  the  course  of  computation,  this  card 
must  be  first.  The  following  cards  contain  two  arrays  of  time  (TB) 
in  sec  and  corresponding  values  of  automatic  ballasting  rate  (BB)  in 
lbs/sec.  The  rules  for  writing  these  data  as  well  as  their  interpreta¬ 
tion  by  the  program  are  similar  to  those  of  valving  data  (type  8) .  , 

3 

Type  10  data  contain  the  air  density  (RH0O)  in  lbs/ft  at  ground 
level  on  the  first  card.  The  following  cards  contain  two  arrays  of 
altitude  above  sea  levelQ(ZZ)  in  ft  and  the  corresponding  values  of 
air  temperature  (TZ)  in  R.  They  must  be  in  order  of  increasing  ZZ. 

The  number  of  cards  must  not  exceed  100.  ZZ(1)  must  be  equal  to  ground 
level.  ZZ(NT)  must  be  equal  to  or  greater  than  the  expected  maximum 
altitude  of  the  balloon.  Atmospheric  data  corresponding  to  ZZ(NT)  are 
used  whenever  the  balloon  goes  above  ZZ(NT)  (see  subroutine  RHOT) . 


*  See  Section  III-C1  for  an  explanation. 
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Type  11  data  contain  two  arrays  of  altitude  above  sea  level  (ZIR) 
in  ft  and  the  corresponding  values  of  the  equilibrium  radiation  tempera¬ 
ture  of  a  black  ball  (RT1R)  in  °R.  The  rules  for  writing  these  data  are 
the  same  as  those  for  type  10  data.  Whenever  the  balloon  goes  above 
ZIR(NIR) ,  infrareddata  corresponding  to  ZIR(NIR)  are  used  (see  sub¬ 
routine  YPRIME) . 

Type  12  data  contain  two  arrays  of  observed  time  (ETIME)  in  sec  and 
corresponding  values  of  observed  balloon  altitude  above  sea  level  (EALT) 
in  ft.  These  data  are  used  for  flight  correlation.  The  only  restriction 
on  these  data  is  that  they  must  be  in  order  of  increasing  ETIME  and  the 
number  of  cards  is  limited  to  400. 

Finally,  type  13  data  contain  two  arrays  of  time  (TIR)  in  sec  and 
the  corresponding  values  of  the  so  called,  "infrared  factor"  (CIR) . 

These  data  pertain  to  infrared  radiation  from  the  earth  and  the  atmos¬ 
phere.  For  a  predictable  atmosphere,  infrared  radiation  is  computed  by 
the  program  from  the  input  data  (type  11)  of  equilibrium  radiation  tem¬ 
perature  of  a  black  ball.  The  infrared  factor  accounts  for  possible 
transient  effects.  The  rules  of  writing  these  data  are  similar  to  those 
of  "solar  radiation  factor"  data  (type  7). 

For  an  example  of  a  complete  set  of  values  of  input  data  see  Sec¬ 
tion  III-E . 

Multiple  Flights.  More  than  one  flights  can  be  computed  in  one  run. 
This  is  accomplished  by  placing  the  decks  of  input  cards  for  each  flight 
in  consecutive  order.  However,  between  flights  the  control  card  with 
the  number  2  must  be  replaced  by  a  blank  card.  If  the  input  data  of 
certain  type(s)  of  a  following  flight  are  the  sama  as  those  of  the  pre¬ 
ceding  flight,  it  is  not  required  to  include  these  data  in  the  deck  of 
the  following  flight.  It  should  be  emphasized  that  this  holds  for  en¬ 
tire  type(s)  of  data  and  not  for  part  of  the  data  in  a  given  type. 

C.  A  BRIEF  DESCRIPTION  OF  THE  PROGRAM 


The  program  is  composed  of  a  MAIN  routine  and  four  subroutines 
(RNGKTA,  YPRIME,  RHOT  and  PLOT),  which  perform  special  tasks. 

1 .  The  MAIN  Routine 


The  MAIN  reads  the  input  data  from  the  input  tape.  When  a  run 
is  made  with  more  than  one  flights,  the  computer  will  read  the  input 
data  and  compute  one  flight  at  a  time.  The  cards  containing  the  input 
data  as  well  as  the  control  cards  must  be  punched  correctly  and  must  be 
in  the  proper  order,  as  described  in  Section  III-B.  Otherwise,  the  first 
time  that  the  MAIN  finds  an  error,  it  will  terminate  the  entire  run  after 
writing  a  comment  in  the  output  tape,  which  will  help  one  to  locate  the 
error  in  the  data. 
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As  a  further  check  that  the  correct  data  are  used,  the  MAIN 
will  write  in  the  output  tape  all  the  data  used  for  the  flight  before 
proceeding  to  the  computation  of  the  flight.  The  input  data  printed 
in  the  output  should  be  examined  carefully  to  make  sure  that  they  are 
identical  to  the  input  data  fed  in  the  computer. 

A  flight  will  be  terminated  automatically  when: 

(a)  The  entire  payload,  including  the  gondola,  has  been 
dropped.  This  could  happen  by  excessive  "manual" 
ballasting  (see  below)  .  A  comment  will  be  written  ir. 
the  output  tape  indicating  that  this  has  happened. 

(b)  The  balloon  hits  the  ground. 

(c)  The  time  exceeds  the  final  time  of  the  flight  specified 
in  the  input  data. 

Then  the  program  will  proceed  to  the  next  flight. 

A  flight  can  be  terminated  instantly  by  turning  SENSE  SWITCH  6 
on.  The  computer  will  pause.  Turn  this  switch  off  and  press  START 
to  proceed  to  the  next  flight. 

The  MAIN  does  a  considerable  amount  of  computation  involved  in 
the  reduction  of  the  input  data  to  a  suitable  form  and  in  the  required 
interpolations  of  the  input  data.  It,  also,  computes  some  of  the  vari¬ 
ables  in  the  desired  output  form  and  units,  as  shown  in  Section  III-E, 
and  store.;  them  in  the  output  tape.  A  unit  of  output  has  two  lines  of 
print.  The  second  line,  which  contains  the  six  heat  transfer  rates,  is 
obtained  only  when  SENSE  SWITCH  2  is  on. 

The  MAIN  also  stores  internally,  not  in  the  output  tape,  the 
information  needed  by  subroutine  PLOT,  which  is  called  only  once  before 
proceeding  to  the  next  flight.  The  information  stored  consists  of  five 
arrays  which  contain  the  time  (sec)  and  the  corresponding  computed  alti¬ 
tude  (ft)  and  temperatures  (°R)  of  the  fabric  and  gas  as  well  as  the 
temperature  of  air  (°R) .  The  dimension  of  these  arrays  is  400.  For  a 
given  flight,  after  these  arrays  are  fillfed,  no  more  information  is 
stored. 


On-line  printing  of  output  is  controlled  by  SENSE  SWITCH  3. 

With  this  switch  on,  each  time  the  computer  stores  output  in  the  output 
tape  it  also  prints  on-line  the  title  of  the  flight  and  the  following 
information:  Greenwich  Mean  Time,  balloon  altitude  above  sea  level  (ft), 
balloon  velocity  (ft/min),  valving  rate  (in  pounds  of  lost  lift  per 
second),  percent  of  total  gas  valved,  ballasting  rate  (lb/sec),  percent 
of  total  ballast  dropped  and  percent  of  total  gas  exhausted.  Other  in¬ 
formation  is  also  printed  on-line,  as  it  is  required  by  the  program, 
and,  therefore,  the  on-line  print  should  always  be  saved. 
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Automatic  valving  and  ballasting  takes  place  at  rates  specified 
in  the  input  data.  Additional  valving  and  ballasting  ("manual")  can  be 
effected  during  the  course  of  the  computation  by  SENSE  SWITCHES  1  and  4, 
respectively.  When  SENSE  SWITCH  1  is  on,  the  computer  will  valve  gas 
in  addition  to  that  valved  automatically  at  a  rate  (VO)  which  is  speci¬ 
fied  in  the  input  data.  The  same  is  true  with  SENSE  SWITCH  4  and 
ballasting. 

The  MAIN  also  performs  exhausting  of  gas  in  a  rather  artificial 
way.  The  integration  is  carried  out  in  steps  of  time  intervals  equal  to 
H.  At  the  end  of  each  step,  the  instantaneous  volume  of  the  gas,  Y(6), 
is  compared  with  the  inflated  volume  of  the  balloon,  VBM,  as  specified 
by  its  manufacturer.  If  Y(6)  is  equal  to  or  smaller  than  VBM,  no  gas 
is  exhausted  and  the  computer  proceeds  to  the  next  step.  If  Y(6)  is 
greater  than  VBM,  the  exhausting  rate  is  set  equal  to  (Y(6)-VMB)/H  and 
the  same  step  is  repeated.  If  again  Y(6)  comes  out  greater  than  VBM, 
the  exhausting  rate  is  increased  by  (Y(6) -VBM) /H,  and  so  on.  Twenty 
such  iterations  per  integrating  step  are  allowed  automatically.  If  at 
the  end  Y(6)  is  still  greater  than  VBM,  the  computer  will  print  on-line 
"ITERATION  DOES  NOT  CONVERGE,  etc"  and  pause.  In  this  position,  when 
START  is  pressed  with  SENSE  SWITCH  5  on,  the  computer  will  enter  again 
the  iterative  loop  for  more  iterations  and  so  on.  On  the  other  hand, 
when  START  is  pressed  with  SENSE  SWITCH  5  off,  the  computer  will  pro¬ 
ceed  to  the  following  flight.  Our  experience  with  this  program  is  that 
an  adequate  amount  of  gas  is  exhausted  with  only  one  iteration  per  in¬ 
tegrating  step. 

Finally,  the  MAIN  does  an  energy  check,  which  will  be  described 
in  Section  III-D. 

2 .  Subroutine  RNGKTA 


The  actual  integration  of  the  six  differential  equations  is 
carried  out  by  this  subroutine  with  Gill's  fourth  order  Runge  and  Kutta 
method,  which  is  a  stable,  self  starting,  accurate  numerical  integration 
technique.  RNGKTA  requires  another  subroutine,  YPRIME,  in  which  the 
differential  equations  are  stated. 

3.  Subroutine  YPRIME 


YPRIME  contains  the  six  first  order  differential  equations  and 
evaluates  the  first  derivatives  of  the  dependent  variables.  It,  also, 
evaluates  the  various  heat  transfer  rates.  Thus,  YPRIME  is  called  fre¬ 
quently  by  MAIN,  as  well  as  by  RNGKTA,  whenever  information  on  the  de¬ 
rivatives  of  the  dependent  variables  or  the  heat  transfer  rates  is  needed. 

For  the  rate  of  absorption  of  solar  radiation,  YPRIME  computes 
the  intensity  of  solar  radiation  and,  when  sunrise  or  sunset  occurs  at 
the  balloon,  it  computes  the  Greenwich  Mean  Time  and  writes  directly  on 
the  output  tape  the  appropriate  comment  (for  instance,  GMT  =  2.20.13 
SUNSET  AT  BALLOON) . 
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When  the  balloon  goes  above  the  infrared  radiation  field  speci¬ 
fied  in  the  input  data  (type  11),  YPRIME  will  print  on-line  an  approp¬ 
riate  comment  and  the  computer  will  pause.  In  this  position,  if  START 
is  pressed,  the  flight  will  be  continued  using  the  infrared  data  of  the 
highest  altitude  for  as  long  as  the  balloon  remains  above  the  specified 
infrared  field.  On  the  other  hand,  if  START  is  pressed  with  SENSE 
SWITCH  6  on,  the  flight  will  be  terminated. 

In  computing  the  convective  heat  transfer  rates  due  to  the 
atmosphere,  YPRIME  gets  information  about  the  atmosphere  from  subroutine 
RHOT. 


4 .  Subroutine  RHOT 


This  subroutine  computes  the  atmospheric  temperature  and  the 
product  of  atmospheric  temperature  and  density  for  given  altitude.  The 
computation  is  done  according  to  formulae  developed  in  Section  II-F. 

When  the  balloon  goes  above  the  atmospheric  temperature  profile 
specified  in  the  input  data  (type  10),  RHOT  will  print  on-line  an  approp¬ 
riate  comment  and  the  computer  will  pause.  In  this  position,  if  START 
is  pressed,  the  flight  will  be  continued  using  the  atmospheric  data  of 
the  highest  altitude  for  as  long  as  the  balloon  remains  above  the  speci¬ 
fied  atmosphere.  On  the  other  hand,  if  START  is  pressed  with  SENSE 
SWITCH  6  on,  the  flight  will  be  terminated. 

Subroutine  RHOT  is  called  frequently  by  MAIN,  as  well  as  by 
YPRIME,  to  give  information  about  the  atmosphere. 

5 .  Subroutine  PLOT 


After  the  completion  of  the  computation  of  a  flight  and  before 
proceeding  to  the  next  flight,  this  subroutine  is  called  to  store  in 
the  output  tape  information  from  which  the  IBM  1401  will  produce  a  plot 
in  graphical  form.  Subroutine  PLOT  is  furnished  to  the  user  in  binary 
card  form  as  the  FORTRAN  source  deck  is  not  available  for  all  the  in¬ 
cluded  subroutines. 

The  grid  of  the  graph  is  marked  with  small  crosses.  The 
abscissa  represents  time  elapsed  from  the  beginning  of  the  flight  in 
hours,  and  the  ordinate  represents  altitude  above  sea  level  in  feet 
and  temperature  in  °R.  There  are  two  compilations  of  this  subroutine: 
in  one,  the  range  of  the  ordinate  is  from  0  to  147,500  feet  for  altitude 
and  from  300  to  595°R  for  temperature,  and  in  another  from  0  to  100,000 
feet  for  altitude  and  from  300  to  550°R  for  temperature. 

On  the  produced  graph  paper,  points  for  the  computed  and 
observed  balloon  altitude  and  the  computed  temperatures  of  the  air, 
fabric  and  gas  are  marked  with  symbols  described  in  the  following  key: 
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DEWA 


(68) 


C  T  w 

p  ab  b 


where  is  the  specific  heat  at  constant  pressure  of  the  gas. 

Notice  that  Equation  63  relates  the  rates  at  which  the  various  forms 
of  energy  of  the  gas-fabric  system  are  exchanged.  The  interpretation  of 
the  various  terms  is  as  follows.  ETOT  is  the  internal  energy.  DKE  and 
DPE  are  the  rates  of  change  of  kinetic  and  potential  energy,  respectively. 
DDRAG  is  the  rate  at  which  energy  is  expended  to  overcome  drag.  DEWA  is 
the  rate  of  work  done  (on  the  atmosphere)  when  the  balloon  expands.  DEVE 
is  the  rate  of  energy  lost  when  gas  is  valved  and/or  exhausted.  Finally, 
the  sum  of  the  q  terms  is  the  rate  of  net  heat  received  by  the  system. 

Let  the  right  hand  side  of  Equation  63  be  denoted  by  ER.  Integrating 
this  equation  from  time  0  to  t,  one  obtains: 

t 

ENET  *  ETOT  -  ETOTl  *  f  ERdt  (70) 


where  ETOTl  is  the  value  of  ETOT  at  zero  time.  Let  the  integral  of  Equa¬ 
tion  70  be  approximated  by  the  following  first  order  formula: 


where  ER^  is  the  value  of  ER  at  the  end  of  each  time  step  H. 

At  the  end  of  each  Integration  by  RNGKTA,  the  MAIN  routine  evaluates 
ETOT,  and,  therefore,  ENET  from  the  integrated  variables.  Similarly,  it 
evaluates  the  increment  In  EP,  which  is  H(ER.) ,  due  to  each  step  and  keeps 
track  of  the  total  EP.  Clearly,  ENET  and  EP,  thus  evaluated,  are  not  equal 
for  two  reasons  (if  all  integrations  were  exact,  ENET  and  EP  would  be 
identical) .  Firstly,  there  is  an  error  involved  in  the  integrations  done 
by  RNGKTA.  Secondly,  there  is  another  error  involved  in  the  above  approx¬ 
imate  first  order  integration  of  ER.  To  be  sure  the  second  error  is  larger 
than  the  first  one.  Nevertheless,  the  ratio: 

ENET-EP  /ton 

ECBK '  w-  C72) 


which  is  evaluated  by  MAIN  and  stored  in  the  output  tape  each  time  output 
is  stored  there,  has  some  bearing  on  the  accuracy  of  RNGKTA.  If  ECHK  is 
a  small  number,  this  means,  at  least,  that  RNGKTA  is  very  accurate. 
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A  -  Air  Temperature 
B  -  Balloon  Gas  Temperature 
E  -  Observed  Balloon  Altitude 
F  -  Balloon  Fabric  Temperature 
L  -  Overlap  of  A  and  B 
P  -  Overlap  of  A  and  F 
Q  -  Overlap  of  B  and  F 
Z  -  Computed  Balloon  Altitude 

$  -  Overlap  of  E  and  Z,  or  of  A,  B,  E,  F,  Z  and  Grid. 

Since  the  computer  can  mark  a  maximum  of  sixteen  equally  spaced 
points  per  half  hour  of  abscissa,  the  uncertainty  in  time  for  each  marked 
point  can  be  as  large  as  1  minute  and  53  seconds.  Similarly,  there  is  an 
uncertainty  in  the  ordinate.  For  instance,  wi  ■  the  0  to  100,000  ft  com¬ 
pilation  the  computer  can  mark  a  maximum  of  six  points  per  10,000  ft  of 
ordinate  and,  therefore,  the  uncertainty  for  each  marked  altitude  point 
can  be  as  large  as  1,667  ft. 

D.  ACCURACY  OF  INTEGRATION 


With  a  program  as  massive  as  this  program  it  is  highly  desireable 
to  have  a  scheme  by  which  the  accuracy  of  integration  can  be  estimated 
for  given  integrating  time  interval.  With  such  a  scheme  the  integrating 
time  interval  can  be  maximized  for  a  given  desired  accuracy  and,  thus, 
the  total  computer  time  per  flight  can  be  minimi  ed. 

To  this  end,  an  energy  check  is  most  effective.  If  Equation  58 
times  y  and  Equations  59  and  60  are  added  and  the  resulting  equation  is 
rearranged,  the  following  energy  equation  is  obtained: 


—■  (ET0T)  =  -  DKE  +  DPE  -  DDRAG  -  DEWA  +  DEVE  +  92  "  q3  +  q4  +  q5  +  q7  (63) 
where: 


ETOT 

DKE 

DPE 

DDRAG 


Cv  “b  TB  +  CF  "F  TF 


(w_  +  w_  +  w_  +  C_  P .  V_)  — *—  y 
N  G  F  B  B  AB'g-7 


<  PA  VB  ‘  “G  '  "F  -  V  * 


1.21  _  r 

2g  'k  CD 


y2  vb2/3 


(64) 

(65) 

(66) 

(67) 
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where  is  the  specific  heat  at  constant  pressure  of  the  gas. 

Notice  that  Equation  63  relates  the  rates  at  which  the  various  forms 
of  energy  of  the  gas-fabric  system  are  exchanged.  The  interpretation  of 
the  various  terms  is  as  follows.  ETOT  is  the  internal  energy.  DKE  and 
DPE  are  the  rates  of  change  of  kinetic  and  potential  energy,  respectively. 
DDRAG  is  the  rate  at  which  energy  is  expended  to  overcome  drag.  DEWA  is 
the  rate  of  work  done  (on  the  atmosphere)  when  the  balloon  expands.  DEVE 
is  the  rate  of  energy  lost  when  gas  is  valved  and/or  exhausted.  Finally, 
the  sum  of  the  q  terms  is  the  rate  of  net  heat  received  by  the  system. 

Let  the  right  hand  side  of  Equation  63  be  denoted  by  ER.  Integrating 
this  equation  from  time  0  to  t,  one  obtains: 

ENET  *  ETOT  -  ET0T1  ■  f  ERdt  (70) 


where  ET0T1  is  the  value  of  ETOT  at  zero  time.  Let  the  integral  of  Equa¬ 
tion  70  be  approximated  by  the  following  first  order  formula: 


where  ER^  is  the  value  of  ER  at  the  end  of  each  time  step  H. 

At  the  end  of  each  integration  by  RNGKTA,  the  MAIN  routine  evaluates 
ETOT,  and,  therefore,  ENET  from  the  integrated  variables.  Similarly,  it 
evaluates  the  increment  in  EP,  which  is  H(ER.) ,  due  to  each  step  and  keeps 
track  of  the  total  EP.  Clearly,  ENET  and  EP,  thus  evaluated,  are  not  equal 
for  two  reasons  (if  all  integrations  were  exact,  ENET  and  EP  would  be 
identical).  Firstly,  there  is  an  error  involved  in  the  integrations  done 
by  RNGKTA.  Secondly,  there  is  another  error  involved  in  the  above  approx¬ 
imate  first  order  integration  of  ER.  To  be  sure  the  second  error  is  larger 
than  the  first  one.  Nevertheless,  the  ratio: 

•  iHr2  (7?) 


which  is  evaluated  by  MAIN  and  stored  in  the  output  tape  each  time  output 
is  stored  there,  has  some  bearing  on  the  accuracy  of  RNGKTA.  If  ECHK  is 
a  small  number,  this  means,  at  least,  that  RNGKTA  is  very  accurate. 
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Many  flights  have  been  computed  with  H  equal  to  10  and  20  sec.  The 
corresponding  values  of  ECHK  are  of  the  order  of  0.003  and  0.007. 

E.  OUTPUTS 


As  stated  in  Section  III-C,  the  first  part  of  output  is  a  printout 
of  the  input  data  used  for  the  flight.  They  are  printed  in  the  form 
shown  in  Figure  2.  Notice  that  the  appropriate  control  card  identifi¬ 
cation  (IT)  precedes  each  type  of  data.  The  meaning  and  units  of  these 
data  are  as  described  in  Section  III-B.  This  printout  of  data  should 
be  checked  carefully  to  make  sure  that  the  correct  input  data  have  been 
used  by  the  computer. 

The  printout  of  input  data  is  followed  by  the  computed  output  in 
the  form  shown  in  Figure  3.  A  unit  output  contains  two  lines.  The 
contents  of  the  first  line  are  the  values  of  the  variables  appearing  at 
the  head  of  each  page  of  output.  These  variables  are: 

Greenwich  Mean  Time  (GMT)  in  hours,  minutes  and  seconds. 

Time  from  beginning  of  flight  (TIME)  in  sec. 

Altitude  of  balloon  from  sea  level  (ALTITUDE)  in  ft. 

Velocity  of  balloon  (VEL)  in  ft/min. 

Atmospheric  pressure  (PRESS)  in  mbar. 

Atmospheric  temperature  (TA)  in  °R. 

Balloon  fabric  temperature  (TF)  in  R. 

Balloon  gas  temperature  (TB)  in  °R. 

Balloon  volume  (VOLUME)  in  ft  . 

Weight  of  payload  and  balloon  fabric  (LOAD)  in  lb. 

Weight  of  balloon  gas  (GAS  WT)  in  lb. 

Free  lift  (FR  LIFT)  in  lb. 

Cumulative  ballast  dropped  (PERB)  in  percent  of  initial  load 
(payload  plus  balloon  fabric). 

Cumulative  gas  valved  (PERV)  in  percent  of  initial  gas  weight. 

Cumulative  gas  exhausted  (PERE)  in  percent  of  initial  gas  weight. 

Number  of  iterations  (IRS)  required  to  exhaust  an  adequate  amount 
of  gas  during  the  most  recent  integrating  step. 

Energy  check  (ECHK)  as  explained  in  Section  IIL-D. 

The  second  line  of  unit  output  contains  the  six  heating  rates  in 
ftlb/sec.  This  line  appears  in  the  output  only  when  SENSE  SWITCH  2  is 
on.  The  heating  rates  are: 

Rate  of  absorption  of  infrared  radiation  by  balloon  fabric  (IRAB). 

Rate  of  emission  of  infrared  radiation  by  balloon  fabric  (IREM) . 

Forced  convection  heating  rate  of  balloon  fabric  by  air  (FCAF) . 

Natural  convection  heating  rate  of  balloon  fabric  by  air  (NCAF) . 

Natural  convection  heating  rate  of  balloon  fabric  by  balloon 
gas  (NCGF) . 

Rate  of  absorption  of  solar  radiation  by  balloon  fabric  (SLAB), 

When  sunset  or  sunrise  occurs  at  the  balloon,  the  Greenwich  Mean 
Time  of  the  occurrence  together  with  the  appropriate  comment  are  printed 
in  the  output  (see  Figure  3) . 
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STRATOSCOPE  FLIGHT  S4-2  JULY  23«1965 
INPUi  DATA 


IT-  3 

250.00000 

1.50000  1.50000  1.00000  0.69000  0.69000 

0.15000 

IT-  4 

11120.00000 

3320.00000  5505000.00000  1460.00000  0.30000  0.50000 

IT-  5 

542.00000 

542.00000  1000.00000  0.  19.90000  31.75000 

102.60000 

IT-  6 

20.00000 

900.00000  0.  40080.00000  4920.00000  198.39000 

300.00000 

IT-  T 

0. 

l.OoOOO 

'0000.00000 

1.00000 

it-  a 

1.00000 

0. 

-0. 

349SO.OCOOO 

0.30000 

35480.00000 

0.33000 

35980.00000 

0.35000 

36480.00000 

0.38000 

36980. OOOCO 

'  0.43000 

37480. OOOCO 

0.48000 

37980.00000 

0.58000 

38480. OOOCO 

0.70000 

38980.00000 

0.88000 

39480.00000 

1.20000 

39980.00000 

1.50000 

41000.00000 

1.50000 

<7* 

II 

►- 

1.00000 

0. 

0. 

60.00000 

1.41670 

240.00000 

0. 

2040. 000 CO 

1.41670 

2340.00000 

0. 

39780.00000 

3.00000 

40080.00000 

•0. 

IT- 10 

0.07000 

1000.00000 

542.00000 

3400.00000 

535.00000 

5025.00000 

529.00000 

9700.00000 

513.00000 

19400.00000 

483.00000 

25000.00000 

459.00000 

31800.00000 

430.00000 

35950.00000 

413.00000 

40800.00000 

393.00000 

46500. OOOCO 

375.00000 

51000.00000 

360.00000 

54550.00000 

368.00000 

61600.00000 

380.00000 

FIGURE  2  INPUT  DATA  PRINTED  IN  OUTPUT  FOR 

STRATOSCOPE  FLIGHT  S4-2 
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68550. 00000 
79100.00000 
82700.00000 
100000.00000 

IT-11 

1000.00000 

51000.00000 

100000. oooco 

lf»12 

0. 

780.00000 
•  1380.00000 
1980.00000 
2580.00000 
3180. OOOCO 
3780.00000 
4380.00000 
5100.00000 
6780.00000 
8580.00000 
10300.00000 
12100. OOOCO 
13900.00000 
1^700.00000 
17500.00000 
19300. OOOCO 
21180.00000 
22980.00000 
24780.00000 
26580.00000 
28380.00000 
30180.00000 
31980.00000 
33780.00000 
34980.00000 
35580.00000 
36180. OOOCO 
36780.00000 
37380.00000 
37980. 000 CO 
38580.00000 
39180.00030 
39780.00000 
40080.00000 

IT-13 

0. 

41000.00000 


384.00000 

400.00000 

407.00000 

423.00000 


532.00000 

392.00000 

392.00000 


1000.00000 

16000.00000 

25000.00000 

37000.00000 

52000.00000 

59000.00000 

66500.00000 

73000.00000 

80000.00000 

80000.00000 

80000.00000 

80000.00000 

80000.00000 

79500.00000 

78500.00000 

78000.00000 

78000.00000 

78000.00000 

78000.00000 

78000.00000 

77500.00000 

77500.00000 

77500.00000 

77500.00000 

77000.00000 

77000.00000 

75000.00000 

73500.00000 

71000.00000 

68000.00000 

65000.00000 

60000.00000 

56000.00000 

47000.00000 

41500.00000 


1.00000 

1.00000 


FIGURE  2  cont'd.  INPUT  DATA  PRINTED  IN  OUTPUT  FOR 

STR  VTOSCOPE  FLIGHT  S4-2 
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FIGURE  3  SAMPLE  OF  PRINTED  OUTPUT 


Following  the  printed  output,  the  output  contains  a  plot  as  shown 
in  Figure  4.  The  contents  of  this  plot  have  already  been  explained  in 
Section  III-C5,  where  subroutine  PLOT  was  described. 

When  a  flight  is  aborted  because  of  errors  in  the  data,  the  output 
contains  comments  which  will  help  one  to  locate  the  errors. 

F.  OTHER  OPERATING  INSTRUCTIONS 

In  accordance  with  the  requirements  of  the  FORTRAN  MONITOR  SYSTEM 
(FMS) ,  a  few  starred  cards  must  be  used  in  the  order  indicated  below 
when  the  program  is  run  with  FMS. 

* IDENTIFICATION  CARD 
*XEQ 

MAIN  ROUTINE  DECK 
SUBROUTINE  RNGKTA  DECK 
SUBROUTINE  YPRIME  DECK 
SUBROUTINE  RHOT  DECK 
SUBROUTINE  PLOT  DECK 
*DATA 

INPUT  DATA  DECK 

The  asterisk  of  the  starred  cards  must  be  in  column  1,  and  their  con¬ 
tents  must  start  from  column  7. 

The  tapes  that  are  used  in  FMS  are  A2  (logical  5)  for  input  and 
A3  (logical  6)  for  output. 

The  functions  of  the  six  sense  switches  of  FMS,  when  they  are  on, 
have  been  described  in  Section  III-C.  In  summary,  the  functions  of 
these  switches  are: 

SSI:  "Manual"  valving  of  gas,  in  addition  to  automatic 
valving,  at  rate  VO. 

SS2:  Heating  rates  are  printed  in  output. 

SS3:  On-line  print  of  output. 

SS4:  "Manual"  ballasting,  in  addition  to  automatic  ballasting, 
at  rate  BO. 

SS5:  Iterations  for  exhausting  gas  will  be  continued. 

SS6:  Present  flight  is  discontinued. 
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IV. 


CORRELATION  OF  COMPUTED  AND  ACTUAL  FL IGHTS 


A.  A  BRIEF  DISCUSSION  OF  THE  APPROXIMATIONS  IN  THE  MODEL  AND  OF  THE 
UNCERTAINTIES  IN  THE  INPUT  DATA 


To  appreciate  the  correlation  of  a  computed  with  an  actual  flight, 
one  should  have  a  good  understanding  of  the  mathematical  model  and  of 
the  computer  program.  Then  one  would  be  in  a  position  to  know  the  de¬ 
gree  of  correlation  he  should  expect.  In  developing  the  mathematical 
model  we  have  made  many  simplifications,  come  of  them  quite  gross,  of 
the  actual  balloon.  In  addition,  the  computer  program  requires  a  certain 
amount  of  specified  quantities  (input  data),  and  the  question  arises  as 
to  how  accurately  these  data  should  and  could  be  specified.  In  this 
section,  we  will  review  briefly  the  approximations  contained  in  the 
mathematical  model  and  the  factors  which,  according  to  the  model,  play 
a  very  important  role  in  balloon  performance.  Then,  in  the  following 
section,  we  will  present  the  correlation  of  some  recent  balloon  flights 
with  those  computed  by  our  program. 

In  calculating  toe  dynamic  parameters  of  the  balloon  (drag  and 
apparent  additional  mass)  and  all  the  heat  transfer  rates  (convective 
and  radiant) ,  we  have  considered  the  balloon  to  be  a  sphere  witn  the 
same  volume  as  the  actual  balloon.  To  compensate  for  this  rough 
approximation,  we  have  introduced  correction  factors  whose  values  have 
been  determined  by  correlation  of  actual  balloon  flights  with  an 
acceptable  degree  of  accuracy. 

Thus,  the  drag  coefficient  of  the  balloon,  based  on  the  cross- 
sectional  area  of  the  equivalent  sphere,  has  been  taken  equal  to  0.3, 

50%  smaller  than  the  drag  coefficient  of  a  sphere,  since  it  is  expected 
that  a  balloon  must  have  less  drag  than  the  equivalent  sphere.  The 
coefficient  of  apparent  additional  mass,  based  on  the  mass  of  the  dis¬ 
placed  air,  has  been  taken  equal  to  0.5,  which  is  equal  to  that  of  a 
sphere  in  potential  flow  (viscous  effects  have  been  neglected) . 

Evidently,  this  is  a  reasonable  approximation,  since  this  coefficient 
depends  on  volume  more  than  it  depends  on  cross-section. 

Besides  the  spherical  model  approximation,  more  approximations 
have  been  introduced  in  the  calculation  of  the  convective  heat  transfer 
rates.  For  forced  convection,  it  has  been  assumed  that  a  known  formula 
for  spheres  in  laminar  flow  is  applicable  to  balloons,  even  through  the 
Reynolds  number  of  balloons  is  usually  higher  than  the  Reynolds  numbers 
covered  by  this  formula.  In  fact,  due  to  high  Reynolds  numbeis,  the 
flow  of  air  around  a  balloon  is  quite  turbulent,  which  means  that  the 
heat  transfer  is  greater  than  that  given  by  this  formula.  Thus,  we 
have  taken  the  correction  constant  equal  to  1.5.  For  free  convection, 
it  has  been  assumed  that  known  formulae  which  hold  for  vertical  plates 
and  cylinders  are  applicable  to  balloons  with  the  diameter  of  the 
spherical  balloon  equal  to  the  height  of  the  plate.  The  correction 
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constants  have  been  taken  equal  to  1.5  and  1.0  for  the  air-side  and  gas- 
side,  respectively,  since  it  is  expected  tha*-  the  air-side  must  have 
greater  heat  transfer  than  the  gas-side. 

For  the  radiant  heat  transfers,  we  f  alieve  that  the  spherical  model 
is  quite  accurate  and,  therefore,  we  have  not  introduced  any  correction 
cor.^,  '.nts  explicitly.  However,  corrections  can  be  introduced  through  the 
effective  absorptivities  and  emissivities . 


another  basic  approximation  introduced  in  the  theory  is  that  both 
the  gas  and  the  fabric  are  characterized  by  an  average  temperature.  The  < 
main  support  for  the  validity  of  this  approximation  stems  from  the  fact  ' 

that  balloons  experience  very  violent  rotations,  which  tend  to  mix  the 
gas  as  veil  as  to  equa’ ize  the  temperature  of  the  fabric  and  the  amount 
of  solar  radiation  incident  on  the  fabric  as  a  whole.  However,  any  heat 
transfer  associated  with  the  rotations  of  balloons  has  not  been  considered! 
in  the  mathematical  model.  Heat  transfer  as  well  as  other  effects  * 

associated  with  the  horizontal  drift  of  balloons  have  also  been  neglected.  ; 


As  pointed  out  earlier,  the  exhausting  process  of  the  model  is  quite 
artificial.  This  may  give,  sometimes,  a  bad  correlation  as  the  balloon 
arrives  at  its  ceiling. 


There  &..A  quite  a  few  parameters  which  play  an  important  role  in 
balloon  performance  and  which  must  be  specified  in  the  computer  program 
as  input  data.  One  of  the  most  important  ones  is  the  initial  free  lift. 
It  has  a  pronounced  effect  on  the  rise  of  a  balloon  to  its  ceiling.  The 
initial  free  lift  is  computed  by  subtracting  the  weight  of  the  payload 
from  the  gross  lift  (minus  weight  of  fabric)  which  is  measured  on  ground 
with  a  scale.  Usually,  the  initial  free  lift  is  not  more  than  about  10% 
of  the  gross  lift.  If  there  is  a  1%  error  in  both  the  measurement  of  the 
payload  and  the  measurement  of  the  gross  lift,  the  uncertainty  in  free 
lift  can  be  as  large  as  20%. 


The  initial  values  of  temperature  of  the  balloon  gas  and  fabric  must 
be  specif iea  as  input  data,  ’'or  a  balloon  which  has  lust  been  filled 
with  gas  and  is  waiting  to  be  released,  specifying  these  temperatures  is 
net  an  easy  task.  In  our  practice  of  the  computer  program  we  have  been 
taking  both  of  these  temperatures  approximately  equal  to  the  air  tempera¬ 
ture  at  ground  level.  These  temperatures  have  an  important  effect  on 
balloon  performance  only  during  the  initial  part  of  the  ascent  phase. 

The  effective  absorptivity  of  the  balloon  for  solar  radiation  and 
its  effective  emissivity  and  absorptivity  in  infrared  are  three  more 
crucial  factors  which  are  sometimes  difficult  to  specify  very  accurately. 
This  is  so  because  these  parameters  are  effective  for  multiple  passes  of 
ra_iation,  the  balloon  gas  being  considered  transparent.  The  one-pass 
absorptivity  or  emissivity  of  a  fabric  can  be  measured  very  accurately. 
But,  in  order  to  calculate  the  effective  multiple-pass  absorptivity  or 
emissivity,  the  one-pass  reflectivity  of  the  fabric  must  also  be  known 
(see  Section  II.H.l).  Accurate  measurement  of  reflectivity  is  not  an 
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easy  task. 

The  black  ball  equilibrium  radiation  temperature  must  also  be 
specified  accurately,  since  a  2,5%  uncertainty  in  this  temperature  (say, 
10°R  in  a  nominal  temperature  of  400°R)  would  give  an  uncertainty  of  10% 
in  radiation  Intensity.  Infrared  radiation  plays  an  Important  >e  in 
balloon  performance,  especially  in  the  absence  of  solar  radiati  .  as  is 
the  case  quite  often  for  a  major  part  of  the  flight  of  large  balloons. 

In  order  to  specify  this  temperature  accurately,  a  black  ball  flight 
should  be  made  before  launching  a  balloon.  If  such  a  flight  is  not 
possible  and  there  is  no  other  way  of  specifying  this  radiation  field,  as 
for  instance  from  black  ball  measurements  made  at  the  launch  site  or  an¬ 
other  similar  site  under  similar  seasonal  conditions,  the  rule  prescribed 
in  Section  II. H» 2  should  be  followed  in  specifying  this  field.  In  any 
case,  a  change  in  the  radiation  temperature  of  the  above  order  can  easily 
take  place  overnight. 

Air  temperature  data  play  a  very  important  role,  but  usually  they 
can  be  specified  accurately.  The  same  is  true  for  valving  and  ballasting 
data. 


Finally,  it  should  be  pointed  out  that  observed  balloon  altitude  is 
computed  usually  from  pressure  measurements  using  the  hydrostatic  law  and  a 
standard  atmosphere  for  atmospheric  temperature.  It  can  be  shown  that 
possible  deviations  of  the  atmospheric  temperature  from  that  of  a 
standard  atmosphere  can  give  changes  in  altitude  of  the  order  of  a  couple 
of  thousand  feet  for  altitudes  above  25,000  ft.  or  so. 

B .  FLIGHT  CORRELATION 


Three  recent  flights  are  correlated  and  discussed  in  this  section. 
Other  flights  have  been  treated  in  our  previous  report  (Ref,  2).  These 
three  flights  include  a  variety  of  effects:  sunset,  sunrise,  ascent, 
overnight  performance,  descent,  valving,  ballasting,  different  fabric 
materials,  etc.  Fairly  accurate  and  reliable  input  data  can  be  prescribed 
for  all  of  them.  Keeping  in  mind  the  discussion  in  the  preceding  section, 
one  can  say  that  the  correlation  of  computed  and  observed  balloon  tra¬ 
jectories  is  reasonable  and  in  fact  very  good. 

The  correlation  of  these  and  other  flights  was  obtained  with  the 
following  values  of  the  five  correction  constants: 


CB  -  0.5 
CD  -  0.3 
C4  *  1.5 
C5  -  1.5 
C6  -  1.0 


These  values  are  compatible  with  theoretical  considerations  as  it  was 
pointed  out  in  the  preceding  se  tion.  We  recommend  that  these  values  be 
used  in  predicting  balloon  performance. 
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1 .  Stargazer  Manned  Flight 


This  balloon  was  launched  from  Holloman  Air  Force  Base,  New 
Mexico  on  December  13,  1962.  It  was  launched  at  1100  MST,  reached  its 
ceiling  of  81,000  ft.  in  about  90  minutes  and  dropped  to  71,000  ft.  after 
sunset  which  occurred  about  6  hours  and  20  minutes  after  launch. 

The  input  data  used  for  this  flight  are  shown  in  Figure  5.  The 
fabric  of  the  balloon  was  GT-12  (mylar  and  scrim),  which  has  a  specific 
heat  of  250  ft  lb  /lb  °R.  Notice  the  value  of  0,69  for  ABIR  and  EMIR 
and  the  value  of  0.15  for  ABUV.  These  values  were  estimated  by  our 
computer  program  (see  Section  II.H.l)  and  then  were  corrected  somewhat 
for  best  flight  correlation.  Notice,  also  that  the  initial  gas  and 
fabric  temperatures  (519°R)  are  9°R  above  the  air  temperature  at  ground 
level.  The  initial  free  lift  specified  in  the  flight  data  was  6%  of 
the  gross  load.  To  obtain  a  reasonable  correlation  in  the  ascent  phase, 
the  initial  free  lift  has  been  increased  to  9.6%  of  gross  load  or  650  lb. 
The  infrared  radiation  field  (IT  =  11)  is  specified  according  to  the 
approximate  rule  given  in  Section  II. H. 2.  Ballasting  and  valving  are 
identical  to  those  of  the  actual  flight.  The  air  temperature  data  were 
given  by  the  Holloman  base. 

The  correlation  of  computed  and  actual  flights  is  shown  in 
Figure  6.  Notice  that  the  agreement  is  excellent  not  only  in  the 
ascent  phase  but  also  during  and  after  sunset.  In  the  computed  flight, 
14.5%  of  the  initial  gas  was  exhausted  by  th?.  time  the  balloon  was 
stabilized  at  its  ceiling.  An  additional  0.7%  was  exhausted  when,  in 
anticipation  of  sunset,  ballasting  took  place  a  little  prematurely.  Thus, 
a  total  of  15.2%  of  gas  was  exhausted  as  compared  to  0.9%  valved  as  the 
balloon  was  reaching  its  ceiling.  The  amount  of  ballast  dropped  during 
sunset  was  7.2 %  of  gross  load. 

2 .  Thermistor  Flight 

This  flight  was  sponsored  by  ONR  and  was  launched  by  NCAR  from 
Page,  Arizona,  on  October  18,  1964.  Arthur  D.  Little,  Inc.,  requested 
this  flight  in  order  to  measure  the  temperature  of  the  balloon  gas  as 
well  as  the  temperature  of  air  at  various  distances  from  the  balloon. 

This  was  accomplished  by  various  thermistors  which  were  placed  by 
Arthur  D.  Little,  Inc.,  in  and  near  the  balloon.  The  principal  objective 
of  this  flight  was  to  acquire  knowledge  about  the  temperature  of  the 
balloon  gas  and  to  compare  this  temperature  with  gas  temperatures  computed 
by  our  program.  Such  a  comparison  is  the  ultimate  check  of  the  validity 
and  accuracy  of  our  model.  A  description  of  the  flight  and  the  instrument¬ 
ation  will  be  published  in  the  near  future. 

A  print-out  of  the  input,  data  is  shown  in  Figure  7.  The  fabric 
of  the  balloon  was  1.5  mil  polyethelene  film.  Its  specific  heat  is  equal 
to  428  ft  lb  /lb  °R.  Ths  radiative  parameters  of  this  fabric  (0.17  and 
0,10)  have  been  estimated  with  our  computer  program  (see  Section  II.H.l). 
The  initial  temperature  of  the  balloon  gas  (519°H)  and  of  the  balloon 
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8H  1  STARGAZER  MANNEO  FLIGHT 
INPUT  DATA 


12-13-63 


1800  GMT 


IT-  3 

250. OOOCO 

1.50000 

1.50000 

1 .00000 

0.69000 

0.69000 

0.15000 

IT-  4 

4800. OOOCO 

1966.00000 

3120000.00000 

650.00000 

0.30000 

0.50000 

IT-  5 

519. OOOCO 

519.00000 

4290.00000 

0. 

-23.19000 

32.91000 

105.94000 

IT-  6 

10.00000 

900.00000 

0. 

43200.00000 

64800.00000 

91.36000 

300.00000 

IT-  7 

-0. 

50000. OOC CO 

1.00000 
l. 00000 

IT-  8 

1.00000 

0. 

4500.00000 
4607.00000 
6000.00000 
6111.00000 
44700.00000 
44916.00000 
45600.00000 
45816. OOOCO 
48300.00000 
48517.00000 
50400.00000 
50490. OOOCO 
51300.00000 
51420.00000 
52800.00000 
52854.00000 
53640.00000 
53698.00000 
57600.000  CO 
57660.00000 
59280.00000 
59336.00000 
59820.00000 
59875.00000 
80000.00000 

IT-  9 

1.00000 

0.  0. 

19260.00000  1.00000 

19335. OOOCO  0. 

20340.00000  1.00000 

20380.00000  0. 

21600.00000  1.00000 

21675.00000  0. 

22560.00000  1.00000 

22635.00000  0. 

22740. OOOCO  1.00000 

22820.00000  0. 


23460.00000  l. 00000 

23540.00000  0. 

24120.00000  1.00000 

24160. 000^0  0. 

24600.00000  1.00000 

24640.00000  0. 

61500.00000  1.00000 

61580. OOOCO  0. 


FIGURE  5  INPUT  DATA  FOR  STARGAZER  MANNED  FLIGHT 
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61680.00000 
61955.00000 
61800.00000 
61 963.00000 
61980.00000 
61955.00000 
62280.00000 
62320. 000C0 
62340.00000 
62495.00000 
63180.00000 
64220.00000 
65100. OOOCO 
65105.00000 
65760.00000 
65765.00000 
1 00000. OOCOO 

IT*  10 

0.06694 
4290.00000 
6571.39956 
39498.00000 
52696.00000 
56000.00000 
59304. OOOCO 
72529.00000 
105660.00000 

IT-11 

4290.00000 
39498. OOOCO 
100000.00000 

IT-12 

200. OOOCO 
1400.00000 
2600.00000 
3800.00000 
5000.00000 
6200. OOOCO 
7400.00000 
8600.00000 
9800.00000 
11000. OOOCO 
12200.00000 
13400.00000 
14600. OOOCO 
15800.00000 
17000.00000 
18200.00000 
19400. OOOCO 
20600.00000 


21800.00000 
23000.00000 
24200.00000 
25400.00000 
26600. OOOCO 
27800.00000 
29000.00000 
30200. OOOCO 
31400.00000 
32600.00000 

IT-13 

-0. 

50000.00000 


FIGURE  5  cont' 


1.00000 

0. 

1.00000 

0. 

1.00000 

0. 

1.00000 

0. 

1. 00000 

0. 

1.00000 

0. 

1.09000 

0. 

0.50000 

0. 

0. 


509.85000 

506.15000 

389.07000 

370.35000 

365.67000 

365.67000 

383.67000 

419.67000 


500.00000 

392.00000 

392.00000 


5000.00000 
23500.00000 
40000.00000 
57500.00000 
73500.00000 
00500.00000 
0 1000. 00000 
81000.00000 
81000.00000 
81000.00000 
81000.00000 
81000.00000 
81000.00000 
81000.00000 
81000.00000 
81000.00000 
80500.00000 
80400.00000 


80000.00000 

79500.00000 

78000.00000 

76500.00000 

74500.00000 

73000.00000 

72000.00000 

71700.00000 

71300.00000 

71000.00000 


1.00000 

1.00000 
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*♦♦♦♦•♦♦♦«*»  Hi 


+  *«>**♦♦**♦»♦♦*♦ 


'44444*4444A*4*4 


444*44444444444! 


l 


44444*4444444*4 


*************** 

E 


♦♦♦♦*♦♦♦♦*♦♦♦♦♦*< 


mHHHHHl* 

Z 


l  A  F  F 

a  a  r 


i*****4**44A*444 


90000 

929 


*********************  t*tt**t»i**t*tf*  t**i*t*$$4*t**t*$*»t*  ft*  tt**i 

*  t  *  * 


•0000 

900 


70000 

479 


»i«Hii»HH  44  44444*  ***********  *******  ******************  ********* 
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fabric  (525°R)  were  deduced  from  the  measurements  of  some  thermistors. 
Notice  that  the  gas  is  colder  than  the  fabric,  the  latter  being  at 
ground  air  temperature.  The  air  temperature  data  are  actual  temperature 
measurements  made  by  a  thermistor  hanging  below  the  gondola.  The  infrared 
radiation  data  are  specified  according  to  the  approximate  rule  given  in 
Section  II. H. 2.  The  observed  altitude  data  (IT  =  12)  were  computed  from 
measurements  of  both  atmospheric  temperature  and  pressure  using  the 
hydrostatic  law.  Therefore,  we  believe  that  they  are  very  accurate. 

The  correlation  of  computed  and  actual  flights  is  excellent  as 
shown  in  Figure  8.  Unfortunately,  the  thermistor  measuring  air  temper¬ 
ature  ceased  to  function  properly  after  the  balloon  reached  its  ceiling 
and,  therefore,  we  cannot  continue  the  correlation  to  the  end  of  the 
ceiling  phase.  Four  hours  and  twenty  minutes  after  launch,  the  cut- 
down  command  was  given  and  the  flight  was  terminated.  The  payload  was 
brought  to  ground  by  parachute. 

The  correlation  of  computed  and  measured  temperatures  is  shown 
in  Figure  9.  Several  thermistors  were  suspended  inside  the  balloon  along 
its  centerline.  Only  the  upper  three  functioned  properly  throughout  the 
entire  flight.  They  were  at  a  fixed  distance  of  7,  19  and  31  feet  from 
the  top  of  the  balloon.  The  relation  of  the  temperature  of  each 
thermistor  to  the  average  temperatures  of  the  gas  and  fabric  depends, 
of  course,  on  the  position  of  the  thermistor  with  respect  to  the  volume 
of  the  balloon.  The  position  of  each  thermistor  during  the  flight  can 
be  deduced  roughly  from  the  sketch  on  the  lower  right  of  Figure  9,  in 
which  the  balloon  is  represented  grossly  by  a  sphere.  The  diameter  of 
the  balloon  was  about  26,  46,  and  80  feet  at  ground,  tropopause  and 
ceiling,  respectively.  Thus,  at  ground,  the  lower  thermistor  was  buried 
in  the  loose  fabric. 

With  the  aid  of  this  sketch,  a  careful  examination  of  Figure  9 
shows  that  the  correlation  of  computed  average  temperatures  and  measured 
temperatures  is  very  good.  The  maximum  deviation,  at  some  parts  of  the 
flight,  is  less  than  10  R.  During  most  of  the  flight,  the  deviation  is 
much  less  than  1U  R. 

We  think  that  this  flight  has  given  considerable  support  to  the 
validity  of  our  analysis,  and  that  it  has  verified  our  long  held  view 
that  the  fabric  of  balloons  experiences  rather  low  temperatures  through 
tropopause.  The  temperature  of  the  top  thermistor  must  be  closely 
coupled  to  the  fabric  temperature.  As  shown  in  Figure  9,  the  fabric  of 
this  balloon  was  at  a  temperature  of  about  380  R  (-78  F)  at  tropopause. 

The  small  oscilations  in  the  computed  temperatures  at  tropopause 
are  due  to  corresponding  oscillations  in  the  input  air  temperature  data. 

3 .  Stratoscope  Flight  S4-2 

Stratoscope  II  Flight  S-4  (Photo)  was  launched  from  Palestine, 
Texas,  on  July  23,  1965,  at  122  GMT.  It  reached  its  ceiling  of  80,000  ft 
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FIGURE  8  PLOTTED  OUTPUT  FOR  THERMISTOR  FLIGHT 


after  one  hour  and  twenty- five  minutes,  approximately  35  min  after  balloon 
sunset.  It  dropped  about  3,000  ft  during  the  night.  It  was  brought 
down  the  next  day  during  sunrise. 

The  input  data  are  shown  in  Figure  2.  The  fabric  of  the  balloon 
was  GT-12  (mylar  and  scrim) .  Notice  that  the  fabric  properties  are 
identical  to  those  of  the  Starr  zer  Manned  Flight.  The  initial  free 
lift  (1,460  lb)  is  exactly  the  value  specified  in  the  flight  data.  The 
initial  balloon  gas  and  fabric  temperatures  (542  R)  are  taken  equal  to 
air  temperature  at  ground  level.  The  air  temperature  data  are  not 
measured  values  for  Palestine.  They  were  computed  from  the  table  of 
Reference  12  pertaining  to  the  July  subtropical  atmosphere  (30  N) .  The 
infrared  data  have  been  computed  according  to  the  approximate  rule  given 
in  Section  II-H.2. 

The  correlation  of  actual  and  computed  flights  is  shown  in 
Figure  4.  Considering  the  fact  that  standard  instead  of  actual  atmos¬ 
pheric  data  have  been  used,  the  correlation  is  fairly  good  during  the 
ascent  and  early  ceiling  phases.  From  this  correlation  and  from  the 
excellent  correlation  of  the  Stargazer  Manned  Flight,  it  appears  that 
the  following  values  for  the  properties  of  the  GT-12  (mylar  and  scrim) 
fabric  are  acceptable: 

Cl  *  250  ftlb/lb°R 
ABIR  -=  0.69 
EMIR  -  0.69 
ABUV  =  0.15 

A  better  correlation  during  the  latter  part  of  the  ceiling  phase 
can  be  obtained  by  reducing  the  infrared  radiation  during  the  night.  This 
can  be  achieved  through  the  infrared  radiation  factor  data  (IT  *  13). 
Leaving  all  the  other  data  as  they  are  and  changing  the  infrared  radia¬ 
tion  factor  data  as  follows: 

IT  *  13 


0.0 

1.0 

7200.0 

1.0 

36000.0 

0.9 

41000.0 

1.0 

we  obtain  the  correlation  shown  in  Figure  10.  The  above  data  mean  that 
from  2  to  10  hours  after  launch  the  Infrared  radiation  incident  on  the 
balloon  is  reduced  by  ig%.  This  means  that  the  black  ball  radiation 
temperature  is  about  10°R  less  than  that  given  in  the  input  data  (IT  -  11). 
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FIGURE  10  PLOTTED  OUTPUT  FOR  STRATOSCOPE  FLIGHT  S4-2  WITH 
REDUCED  (by  10%)  INFRARED  RADIATION  DURING  NIGHT 
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FIGURE  lOcont’d.  PLOTTED  OUTPUT  FOR  STRATOSCOPE  FLIGHT  S4-2  WITH 
REDUCED  (by  10%)  INFRARED  RADIATION  DURING  NIGHT 


Srtlfur 


Such  a  change  in  the  infrared  radiation  field  from  day  to  night  has,  in 
fact,  been  noticed  (see  Section  II-H.2).  Notice  that  the  correlation 
during  the  night  is  now  much  better.  Evidently,  reduction  of  infrared 
radiation  caused  the  actual  balloon  to  drop  3,000  ft  during  the  night. 

The  correlation  of  the  last  part  of  the  desent  phase  is  not 
good.  This  must  be  due  to  the  fact  that  there  is  some  uncertainty  in 
the  time  of  closing  of  the  helium  valve. 
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LIST  OF  SYMBOLS 


u 


03 

,o 

O 

SO 

< 

Fortran 

Description 

Units 

A 

Effective  area  of  balloon  in  direction  of  motion, 
cross-sectional  area  of  spherical  balloon. 

ft2 

b 

n 

TZB(I) 

See  Equation  25. 

°R 

B 

B 

Ballasting  rate. 

lb/sec 

CB 

CB 

Apparent  additional  mass  coefficient  of  balloon. 

CD 

CD 

Drag  coefficient  of  balloon. 

CF 

Cl 

Specific  heat  of  balloon  fabric. 

ftlb/lb°R 

CP 

CP 

Specific  heat  of  balloon  gas  at  constant  pressure. 

f tlb/lb°R 

cv 

CV 

Specific  heat  of  balloon  gas  at  constant  volume. 

f tlb/lb°R 

C4 

C4 

Correction  constant  for  forced  convection  heat 
transfer  from  air  to  balloon  fabric  . 

C, 

C5 

Correction  constant  for  free  convection  heat 
transfer  from  air  to  balloon  fabric. 

C6 

C6 

Correction  constant  for  free  convection  heat 
transfer  from  balloon  gas  to  balloon  fabric. 

D 

Diameter  of  spherical  balloon. 

ft 

DDR4G 

DDRAG 

Rate  of  energy  spent  to  overcome  drag. 

ftlb/sec 

DEC 

SEC 

Declination  of  sun. 

deg 

DEVE 

DEVfi 

Rate  of  energy  lost  when  gas  is  valved  and/or 
exhausted . 

ftlb/sec 

DEWA 

DEWA 

Rate  of  work  done  by  balloon  gas  on  atmosphere 
when  balloon  expands. 

ftlb/sec 

DKE 

DKE 

Rate  of  change  of  kinetic  energy  of  balloon. 

ftlb/sec 

DPE 

DPE 

Rate  of  change  of  potential  energy  of  balloon. 

ftlb/ sec 

E 

E 

Rate  of  volumetric  gas  flow  due  to  exhausting. 

ft^/sec 

ECHK 

ECHK 

Energy  check  ratio  of  integration  (see  Section  III-D) . 

ENET 

ENET 

Difference  of  internal  energy  of  gas  and  fabric  from 
that  at  beginning  of  flight. 

ftlb 

EP 

EP 

Same  as  ENET,  but  found  from  integration  (see  Section 
III-D) . 

ftlb 

ER 

Rate  of  change  of  mechanical  and  thermal  energy  of 
balloon. 

ftlb/sec 
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o 


« 

u 

J3 

4) 

00 

3 

Fortran 

Description 

Units 

ERi 

Value  of  ER  at  the  end  of  each  step  of  integration. 

ftlb/sec 

ETOT 

ETOT 

Internal  energy  of  balloon  gas  and  fabric. 

ftlb 

ETOT1 

ET0T1 

Value  of  ETOT  at  beginning  of  flight. 

ftlb 

FL 

FL 

Free  lift. 

lb 

g 

G 

Acceleration  of  gravity  («  32.2). 

2 

ft/sec 

GHA 

GHA 

Greenwich  hour  angle  at  beginning  of  flight. 

deg 

GL 

GL 

Gross  lift. 

lb 

h 

Coefficient  of  heat  transfer. 

ftlb/ ft^sec°R 

k 

Thermal  conductivity. 

ftlb/ft3ec°R 

kA 

CAIR 

Thermal  conductivity  of  air. 

ftlb/ftsec°R 

CHE 

Thermal  conductivity  of  balloon  gas. 

ftlb/ftsec°R 

i. 

V 

Lost  lift  per  unit  time  due  to  valving. 

lb/sec 

LAT 

XLAT 

Latitude  of  balloon. 

deg 

LHA 

Local  hour  angle. 

deg 

LONG 

XL0NG 

Longitude  of  balloon. 

deg 

m 

A  TRM 

Air  mass . 

m 

o 

Air  mass  at  sea  level . 

ml 

AIRMl 

Air  mass  for  point  B  in  Figure  1. 

m2 

AIRM2 

Air  mass  corresponding  to  the  supplement  of  actual 

L 

azymuth  angle  of  sun. 

\ 

XMA 

Molecular  weight  of  air. 

XMB 

Mole<.  ar  weight  of  balloon  gas. 

n 

Index. 

P 

P 

Atmospheric  pressure. 

lb/ft2 

Pn 

Atmospheric  pressure  at  altitude  z^. 

lb/ft2 

Po 

PZER 

Atmospheric  pressure  at  ground  level. 

lb/ft2 

pl 

PZ1 

Atmospheric  pressure  at  point  B  in  Figure  1. 

ib/ft2 

PrA 

Prandci  number  of  air. 

Pr 

H 

Prandtl  number  of  balloon  gas. 

DEH2 

Rate  of  absorption  of  infrared  radiation  by  balloon 

ttlb/sec 

i 

fabric. 

DEH3 

Rate  of  emiaaion  of  Infrared  radiation  by  balloon 

ftlb/sec 

J 

fabric. 
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Description 


Units 


« 

U 

X> 

0) 

00 


c 

* 

k- 

u 

u 

o 

Cm 


I 


q4 

DEH4 

Rate  of  heating  of  balloon  fabric  by  air  through 
forced  convection. 

ftlb/sec 

q5 

DEH5 

Rate  of  heating  of  balloon  fabric  by  air  through 
free  convection. 

ftlb/sec 

q6 

DEH6 

Rate  of  heating  of  balloon  fabric  by  balloon  gas 
through  free  convection. 

ftlb/sec 

q7 

DEH7 

Rate  of  absorption  of  solar  radiation  by  balloon 
fabric . 

ftlb/sec 

Q 

Energy  supplied  to  a  system. 

ftlb 

R 

R 

Gas  constant  per  mol  (**  1,545),  or 
average  reflectivity  of  balloon  fabric,  or 

ftlb/lbmol°R 

R 

r 

Re 

RETH 

average  radius  of  earth  («  20,903,520). 

Infrared  average  reflectivity  of  balloon  fabric. 

Reynolds  number. 

ft 

s 

n 

TZA(I) 

See  Equation  24. 

°R/ft 

S 

Surface  .  tea  of  spherical  balloon. 

ft2 

t 

X 

Time  from  beginning  of  flight. 

.  'C 

T 

TRANS 

Intensity  of  solar  radiation  relative  to  that  outside 
atmosphere. 

' 

T 

TA 

Temperature  of  air. 

°«  1  ; 

tb 

Y<4) 

Temperature  of  balloon  gas. 

o_  ■  j 

R  i  : 

tf 

Y(3) 

Temperature  of  balloon  fabric. 

oB  l  1 

■  1 

T 

n 

TZ(I) 

Temperature  of  air  at  altitude  2^. 

;  ! 

T 

r 

Black  ball  equilibrium  radiation  temperature. 

°a  1  1 

U 

Internal  energy  of  a  system. 

ftib  i 

•a  ■  ■  t  | 

• 

V 

Rate  of  volumetric  gas  flow  due  to  valving. 

ft3/s*'  i  1 

VB 

Y<6) 

Volume  of  balloon. 

ft3  ‘  \ 

WB 

Y<5> 

Weight  of  balloon  gas. 

lb  j 

wc 

Weight  of  payload. 

lb 

V 

V 

WF 

Weight  of  balloon  fabric. 

lb 

W 

Work  dona  by  a  system. 

ftlb 

X 

Product  of  Crashof  and  Praodtl  numbers. 
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o 


c 

u 

X 

u 

V 

4-1 

u 

. — i 

o 

< 

Cx* 

Description 


Units 


> 

z 

7 


v(2> , r(i) 
Y(l) 
ZZ(I) 


zx  ZZ1 

4 

4r 

4  ref  ^ 
Of  vRf  ABUV 


<T 

£ 

^ef 


6: 

£ 


r 

ref 


EMIR 


VISA 

>*  H  VISH 

/*A  RH0 

p  B  RH0H 

P  n 

O'  SB0LZ 

'Tv 


Vertical  velocity  of  balloon.  ft/sec 

Altitude  of  balloon  above  sea  level.  ft 

Altitudes  at  which  atmospheric  input  data  are  given.  ft 

Altitude  of  ground  above  sea  level.  ft 

Altitude  of  point  B  in  Figure  1.  ft 

Average  absorptivity  of  balloon  fabric. 

Infrared  average  absorptivity  of  balloon  fabric. 

Effective  infrared  absorptivity  of  balloon  fabric, 
for  multiple  absorption. 


Effective  ultraviolet  absorptivity  of  balloon  fabric 
for  multiple  absorption  (solar  radiation). 

Azymuth  angle  of  sun.  deg 

Average  emissivity  of  balloon  fabric. 

Effective  average  emissivity  of  balloon  fabric  for  ' 
multiple  emission. 

Infrared  average  emissivity  of  balloon  fabric. 

Effective  infrared  emissivity  of  balloon  fabric  for 


multiple  emission. 

Viscosity  of  air.  lb/ftsec 

Viscosity  of  balloon  gas.  lb/ftsec 

3 

Density  of  air.  lb/ft 

3 

Density  of  balloon  gas.  lb/ft 

3 

Density  of  air  at  altitude  z  .  lb/ft 

n  10  £  o 

Stefan  -  Boltzmann  constant  (=  3.6995  x  10  ).  ftlb/ft  sec  R 


Infrared  average  transmissivity  of  balloon  fabric. 
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n 

0 

fe 

Description 

Units 

AH 

Hours  part  of  Greenwich  Mean  Time. 

hours 

AM 

Minutes  part  of  Greenwich  Mean  Time. 

min 

AS 

Seconds  part  of  Greenwich  Mean  Time. 

sec 

BO 

Manual  ballasting  rate. 

lb/sec 

BB(I) 

Array  of  specified  automatic  ballasting  rates 

Ib/sec 

CAM 

Cosine  of  azymuth  angle  of  sun. 

CCIR 

Infrared  radiation  factor. 

CCSI(I) 

Array  of  specified  solar  radiation  factors. 

CLATD 

Cosine  cf  LAT  times  cosine  of  DEC. 

CIS(I) 

Array  of  specified  infrared  radiation  factors. 

CSI 

Solar  radiation  factor. 

DCIRDT(I) 

Array  of  computed  slopes  of  infrared  radiation  factor 

-1 

.  sec 

DIRDZ(I) 

Array  of  computed  slopes  of  FIR(I) . 

ftlb/ft  sec 

DP 

Plotting  time  interval. 

sec 

DPR 

Printing  time  interval. 

sec 

-1 

sec 

DSIDT(I) 

Array  of  computed  slopes  of  solar  radiation  factor. 

EALT(I) 

Array  of  observed  balloon  altitudes  above  sea  level 
corresponding  to  times  ETIME(I) . 

ft 

ETIME(I) 

Array  of  times. 

sec 

FCAF 

Forced  convection  Seating  rate  of  balloon  fabric  by  air.ftlb/sec 

F(2) 

Vertical  acceleration  of  balloon. 

.  i 

ft /sec 

F(3) 

Time  derivative  of  balloon  fabric  temperature  . 

°R/'  sec 

F(4) 

Time  derivative  of  balloon  gas  temperature. 

°R/sec 

F(5) 

Time  derivative  of  balloon  gas  weight. 

lb /sec 

F(6) 

Time  derivative  of  balloon  volume. 

ft3/sec 

FXR(I) 

Array  of  computed  atmospheric  radiation  fluxes. 

ftlb/ft2sec 

FLO 

Ini'  il  free  lift. 

lb 

FLUXIR 

Infrared  radiation  flux. 

ftlb/ft2sec 
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Desci option 


Units 


IRS, LAMBDA 


Value  of  solar  constant  outside  earth's  atmosphere 
(-  96). 

Greenwich  Mean  Time. 

Greenwich  Meat  Time  at  beginning  of  flight. 

Time  interval  for  integration. 

Rate  of  absorption  of  infrared  radiation  by  balloon 
fabric . 

Rate  of  emission  of  infrared  radiation  by  balloon 
fabric. 

Number  of  iterations  in  exhausting  gas. 

Weight  of  payload  and  balloon  fabric . 

Natural  convection  heating  rate  of  balloon  fabric  by 


ftlb/ft  sec 


hr :min;sec 


ftlb/sec 


ftlb/sec 


PM, PRESS 


RTIR(I) 


RIZ(I) 


Natural  convection  heating  rate  of  balloon  fabric  by 
balloon  gas. 

Cumulative  ballast  droped  in  percent  of  initial  load 
(payload  plus  balloon  fabric). 

Cumulative  gas  exhausted  in  percent  of  initial  balloon 
gas  weight. 

Cumulative  gas  valved  in  percent’ of  initial  balloon  gas 
weight . 

One-third  power  of  balloon  volume. 

Atmospheric  pressure. 

Conversion  factor  from  degrees  to  radians  ygjj  J. 

Density  of  air  at  ground  level . 

Product  of  density  and  temperature  of  air 

Array  of  specified  black  ball  equilibrium  radiation 
temperatures. 

Array  of  products  of  density  and  temperature  of  air 
computed  from  atmospheric  data. 

Rate  of  absorption  of  solar  radiation  by  balloon 
fabric . 


ftlb/sec 


rad/deg 

lb/ft3 

lb°R/ft3 


lb  R/ft" 


ftlb/sec 


SLATD 


Sine  of  LAT  times  sine  of  DEC. 
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Description 


Units 


TAA(I) 

Array  of  computed  air  temperatures  stored  for 
subroutine  PLOT, 

°R 

TB(I) 

Array  of  specified  times  for  automatic  ballasting. 

sec 

TBO 

Initial  temperature  of  balloon  gas. 

°R 

TCSI(I) 

Array  of  specified  times  for  solar  radiation  factor. 

sec 

TFO 

Initial  temperature  of  balloon  fabric. 

°R 

TIR(I) 

Array  of  specified  times  for  infrared  radiation 
factor. 

sec 

TV(I) 

Array  of  specified  times  for  automatic  valving. 

sec 

VALD 

Cumulative  weight  of  valved  gas. 

lb 

VBM 

Inflated  volume  of  balloon. 

ft3 

VRL 

Vertical  velocity  of  balloon. 

ft/min 

VO 

Manual  valving  rate  (in  pounds  of  lost  lift  per 
second) . 

lb/sec 

W(I) 

Array  of  specified  automatic  valving  rates  (in 
pounds  of  lost  lift  per  second) . 

lb/sec 

WBO 

Initial  weight  of  balloon  gas. 

lbs 

WFO 

Height  of  balloon  fabric. 

lb 

WG 

Weight  of  payloaa  and  balloon  fabric. 

lb 

WGO 

Initial  weight  of  payload. 

lb 

WT 

Initial  weight  of  payload  and  balloon  fabric. 

lb 

XL  AH 

Local  hour  angle. 

rad 

XP 

Printing  time. 

sec 

XPL 

Plotting  time. 

sec 

XT 

Specified  final  time  of  flight. 

sec 

XO 

Specified  initial  time  of  flight. 

sec 

XXX(I) 

Array  of  times  stored  for  subroutine  PLOT. 

sec 

Y1(I) 

Array  of  computed  balloon  altitudes  stored  for 
subroutine  PLOT. 

ft 

Y20 

Specified  initial  vertical  velocity  of  balloon. 

ft/sec 

Y3(I) 

Array  of  computed  temperatures  of  balloon  fabric 
stored  for  subroutine  PLOT. 

°R 

Arthur 


Description 


Units 


u 

‘tA 

CO 

c 

u 

nS 

£> 

u 

V 

aJ 

00 

U 

t-A 

0 

< 

to 

Y4(I) 


ZIR(I) 

ZO 


Array  of  computed  temperatures  of  balloon  gas  stored 
for  subroutine  PLOT. 


Amy  of  altitudes  above  sea  level  at  which  RTIR(I) 
is  given. 


Specified  initial  altitude  of  balloon  above  sea  level. 


ft 

ft 


( 

f 

I 

i 

i 

| 

| 
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FIGURE  II 


GENERAL  FL( 


bURE  11  GENERAL  FLOW  CHART 


ooooooooo 


ONR-ADL  COMPUTED  VERTICAL  TRAJECTORY  OF  HIGH 
COMPILATION  8H 


SENSE  SWITCH  CONTROLS 
SENSE  SWITCH  I 

SENSE  SWITCH  2 

SENSE  SWITCH  3 

SENSE  SWITCH  4 

SENSE  SWITCH  5 

SENSE  SWITCH  6 

5  FORMAT  (  12  »  3F 10#  5 ) 

10  FORMAT  {  12,  7E10.3) 
15  FORMAT  ( 77H1  YOU  HAVE 


DOWN 

VALVE  GAS  (VO) 
PRINT  DETAIL 
PRINT  ON-LINE 
DROP  BALLAST  (RO) 
CONTINUE  ITERATIONS 
ABORT  FLIGHT 


THROWN  AW*Y  ALL  OF 


ALTITUDE  BALLOONS 
UP 

NORMAL 

NORMAL 

NORMAL 

NORMAL 

STOP  AFTER  20 
NORMAL 


THE  BALLAST  -  -  INCLU 


10ING  THE  GONDOLA. /1H0/1H0) 

20  FORMAT  (  12,  7F10.4) 

25  FORMAT  (23H  DERATION  DOES  NOT  C0NVERGE/54H  SENSE  SWITCH  5  DOWN  A 
IND  START  TO  CONTINUE  I TERAT I0NS/40H  SENSE  SWITCH  5  UP  AND  START  T 
20  END  RUN/lHO/iHOl 

30  FORMAT  (  12,  2F10.4) 

31  FORMAT! I2,2F10#4) 

35  FORMAT  (  7H0  T I ME=F 10. 5 , 4H  Z=F10.2,?H  ZDOT*F 10. 5 , 5H  TF*F10.4, 
15H  TB=F10.4,5H  WB-Fi0.4,5H  VR*E15.8/  7H0  VBM=E15.8t3H  E*E15.8) 

40  FORMAT  (  7H1  T I ME-F 10. 5 , 8H  TV(1)=F10.5> 

45  FORMAT! 15H  INPUT  IS  WRONG/ l HO/ 1H0 ) 

50  FORMAT  ( 3 1 H  CARD  AFTER  TITLE  CARO  IS  WRONG/LHO) 

55  FORMAT ( 7HI  TIME=F10. 5, 8H  TR(l)=Fi0.5S 
60  FORMAT  (9H  NO  TYPE  12, 5H  DATA/lHO ) 

65  FORMAT  (  7H1  T IME*F 10.  5®9H  TCS I ( l ) ~F 10. 5) 

70  FORMAT  (15H  TOO  MANY  TYPE  I2,5H  DATA/lHO) 

75  FORMAT  (20H1  WBOOT  IS  POSITIVE  F10.5) 

80  FORMAT  (20H  IT  AFTER  DATA  TYPE  I2,12H  IS  NEGAT I VE/ IHO ) 

81  FORMAT  »  7HI  T I ME=F 10 . 5 , 8H  TIRI1)*F10.5) 

90  FORMAT  (16H0  TOO  MANY  TYPE  I2.13H  DATA  CARDS.  /IHO  /IHO) 

91  FORMA T ( 5H0GMT -F 3. 0, F3.0,F3.0,10H  Al T I TUDE-E6.0 , l OH  VELOCI T Y*F6.0, 

1 1 OH  VAL  r  AT E=F5 • 2 , 1 3H  PERCENT  VAi -Fb.2* 10H  RAL  RATE»F5.2, 

213H  PERCENT  RAL=F6.?,13H  PERCENT  6XD=F6.2/1H0/1H0) 

993  FORMAT! lhl,5lH  CARD  AFTER  CARO  TYPE  1  OF  PREVIOUS  PROBLEM  IS  WRONG 
I/IHO) 

994  FORMAT  !  IHI,  12A6  ) 

995  FORMAT  (  12A6  ) 

975  FORMAT  !  12,  5F10.5) 

976  FORMAT! IHO, IX, 3HGMT , 6X , 4HT I  ME , 2X , 8HAL T I TUDE * 4X , 3HVEL, 4X,5HPRESS* 3X 
1,2HTA,5X,2HTF,5X,2HTB,5X, 6H VOLUME , 4X,4HL0AD , 3X » AHGAS  WT*2X,7HFR  LI 
2FT,2X,4HPEKB,2X,4HPERV,2X,4HPERE  2X,3HIRS,3X,4HECHK) 

977  FORMAT  ( IHO, 3F 3. 0 , 1 XF7 .0, 2XE7. 0, 2XF 6 .0, 2XF6. 1 ,2XF5. 1 , 2XF5. I  * 2XF 

15.1,2XF8.0,2XF6.0,2XF6.1,2XF7.2,2XF4.1»2XF4. l,2XF4.1 ,2XI2,2XF6.3) 

5001  F0RMAT(4X,5H!RAB*lPEll.4,3X,5HlREM«lPEil.4,3X,5HFCAF«lPEU.4,3X,5H 
INCAF* IPE l 1.4, 3X ,5HNCGF*1PE 1 1.4,3X,5HSLAB*1PE11.4) 

2000  FORMAT! IHO, IIH  INPUT  OATA) 

2001  FORMAT! IHO, 4H  I T* 12  > 

2002  FORMAT! 7F16. 5) 

2005  FORMAT  (2F16.5) 

2006  FORMAT  (F16.5) 

DIMENSION  SCALE ( 60 ) 

DIMENSION  ZZ(IOO),  TZ  C 100)  *  RTZUOOI,  TZA(IOO),  TZBI100I 
^TMENSION  TCSI1100),CCS1(100I,  TVIIOO),  VVIIOO) 


DIMENSION  TB(IOO),  BB(IOO)  ,  XX ( 100 i 
01  MENS  I  ON  YSIOO),  FUQO),  <H100) 

DIMENSION  TITLE  (12) 

DIMENSION  RT  I R  (  100) 

DIMENSION  ZIR(l00),FIR(l00)»DIRCZ(100),0SIDT(lO0) 

DIMENSION  Z UV( 100) ,FUV* 100) tOUVDZ! 100) 

DIMENSION  Y 1  ( 400 )  »  Y3 ( 400 ) »  Y4 ( 400 ) *  XXX ( 400 ) »  TAA ( 400 ) 

DIMENSION  ETIME(400),FALT(400) 

DIMENSION  TIRUOO),CIR(  100)  ,DCIRDT(  100) 

COMMON  X,  Y,  F,  Q 

COMMON  NT,  ll ,  TZ,  RTZ ,  TZA,  TZB 

COMMON  Cl.  C2,  C3,  C4,  C5,  C6,  C7 

COMMON  TA,  RT,  P,  RHO,  PHI,  WG,  OMEGA,  CSI,  E,  V 

COMMON  XMA,  XMB,  G,  CP,  CV,  R,  CD,  VRM  ,WF 

COMMON  DtH  ,D6H3»DFH4  DEH5,DEH6,DEH7' 

COMMON  WFH, DWFH, WFO 
COMMON  CB 

COMMON  Z1R,FIR,DIRDZ,NIR, ZUV , FUV.DUVQZ , NUV, DSIDT 
COMMON  A81R,EMIR* ABuV 
COMMON  M,Yi>hY3,Y4,TAA,XXX 
COMMON  DP,*.ftXO,ETIME»FALT 
COMMON  CCIR 

COMMON  SIATD,ClATD,RAD,GH* , XLONG, PZER ,FSOL , SBOLZ , AH,AM, AS ,GMTS 
FSQL*96« 

SBOLZ  *3 .6995E-10 
XMH=  4. 

XMA=  28.89 
CP*972. 69 
CV*536.73 
R=1545. 

G  *  32.2 

RAD- 3. 14 l 59/ l 80. 

1001  READ  INPUT  TAPE  5,  995 , ( TT TIE (  I ) ,  1*1,12) 

DO  990  MZY*1 ,400 
Y1(MZY)*0. 

Y3 ( MZY )*0. 

Y4(MZY)*0. 

xxx(mzy;*o. 

TAA ( MZY ) *0. 

990  CONTINUE 

PRINT  994,  (TITLE! I ),  1*1,12) 

999  READ  INPUT  TAPE  5,  10,  IT 
IF { IT)  998,998,1000 

998  WRITE  OUTPUT  TAPE  6.994. (TITLE (I). 1*1. 121 
WRITE  OUTPUT  TAPE  6,50 
CALL  EXIT 

1000  GOTO( 100,200,300,400.500,950,700,800,900,600,960,  980, 1400), it 

997  WRITE  OUTPUT  TAPE  6, 994, ( T I TLF ( I) , I » l ,  12) 

WRITE  OUTPUT  TAPE  6,80, IP 
CALL  EXIT 

970  WRITF  OUTPUT  TAPE  6,994, (TITLE ( I )  ,  l *1, 12) 

WRITE  OUTPUT  TAPE  6,60,!P 

CALL  EXIT 

971  WRITE  OUTPUT  TAPE  6 ,994, ( TI TIE  (!) , 1*1* 12 ) 

WRITE  OUTPUT  TAPE  6,70, IP 


CALL  EXIT 

200 

CALL  EXIT 

100 

WRITE 

OUTPUT 

TAPE 

WRITE 

OUTPUT 

TAPE 

1*3 

WRITE 

OUTPUT 

TAPE 

WRITE 

OUTPUT 

TAPE 

1=4 

WRITE 

OUTPUT 

TAPE 

WRITE 

OUTPUT 

TAPE 

1=5 

WRITE 

OUTPUT 

TAPE 

WRITE 

OUTPUT 

TAPE 

1*6 

WRITE 

OUTPUT 

TAPE 

WRITE 

OUTPUT 

TAPE 

1*7 

WRITE 

OUTPUT 

TAPF 

DO  2004  1  =  1,1 

MSI 

2004 

WRITE 

OUTPUT 

TAPE 

1*8 

WRITE 

OUTPUT 

TAPE 

WRITE 

OUTPUT 

TAPE 

DO  2007  1  =  1,1 

MV  I 

2007 

mRITE 

OUTPUT 

TAPE 

1*9 

WRITE 

OUTPUT 

TAPE 

WRITE 

OUTPUT 

TAPE 

DO  2008  1*1,1 

NBI 

2008 

WRITE 

OUTPUT 

TAPE 

1*10 

WRITE 

OUTPUT 

TAPE 

WHITE 

OUTPUT 

TAPE 

DO  2009  1*1, 

NT 

2009 

WRITE 

OUTPUT 

TAPE 

1*11 

WRITE 

OUTPUT 

TAPE 

DO  2010  1*1, 

NIR 

2010 

WRITE 

OUTPUT 

TAPE 

1*12 

WRITE 

OUTPUT 

TAPE 

2011 


2012 

2013 


994, ( TITLE! 11,1=1,12) 

2000 

2001,1 

2002,Cl,C4,C5,C6,ABIR,EPIR,ABUV 

2001,  I 

2002»WG0,WF0»VBM,FL0»CD»CB 

2001,1 

2002,  TB0,TF0,Z0,  Y20,DEC,XLAT,XL0NG 

2001,1 

2OG2,H,DFR,X0,XT,GMTS,GHA,DP 

2001*1 

2005, Tcsim,ccsm ) 

2001,1 

2006,  VO 

2005,  TV! I) ,VV( I ) 

2001,1 

2006,  BO 

2005, TR(n,BRm 

2001,1 

2006,  RHOO 

2005, ZZ ( I ) » TZ (  I ) 

2001,1 

2005, Z I Rl I) »RT  1RI I  ) 

2001,1 


2005,  EH  ME  ( I )  ,E  ALT  ( 1 ) 
2001,1 

2005, T 1R ( I ) ,C IR( n 


DO  2011  1*1 ,NEX 
WRITE  0U1PUT  TAPE 
1-13 

WRITE  OUTPUT  TAPE 
00  2012  I«l,NCIR 
WRITE  OUTPUT  TAPE 
CLATO*CO$F(XLAT»RAD)*COSF!DEO»RAD) 
SLATD-SINFl XLAT*RAD)*SINF<DEC*RAO) 

rm-zzi  i) 

CAU  RHOT  I  Vf  1 )  ,TFMP,RT  I 

PZER«RT»R/XMA 

FL-FLO 


r 

i 


i 


Y(2)*Y20 

mi  =  tfo 

Y  ( 4 )  =  T60 
X  *  xo 
XP=  XO 
XPL»XO 
WF*WFO 
WG*WGO+WFO 
VALDxO 
wr=wG 
LP*  25 
LPC  *1 
I  Y*0 
NPL*0 

CALL  RHOT  (Y(l),  TEMP,  RT) 

RHO  *  RT  /  TEMP 
p=RT'r,/.Y*'* 

TA  *  TfcPP 

Y(5)=(FL*WG)/(XMA»Y(4)/(XMB*TA)-l.  ) 

W80*Y ( 5  ) 

Y(6)*R»Y(5)*Y(4)/(P*XMB) 

CS I =CCS !  ( 1 ) 

CCIR*CIR(  1  ) 

PH  I  =  Y 1 6 )**( 1 . / 3 .  ) 

C  -  INITIALIZING  TOTAL  FNFRGY  - El - 

CALL  YPRIMfc 

DKE*(WG+Y(5)+.5*RH0»Y(6))»Y(2)»F(2)/G 
OPE*  ( RH0*Y(6 )-WG-Y( 5 ) l«Y( 2 ) 

DDRAG*CD*.60375*PHI»PHI*RM0»Y(2)»ABSF(Y(2))»Y(2)/G 

0EWA*P*F<6) 

DEVE*Y(4,»CP*F(5) 

ET0T*CV*Y(5)*Y(41+C1»WF0»Y(3) 

ETOTl*ETOT 

ENET-.O 

EP*0.0 

fcCHK*0. 

101  LAMBDA*0 

IF  (X  -TVCD)  102,103,103 

102  WRITE  OUTPUT  TAPE  6,40,XtTV(l) 

WRITE  OUTPUT  TAPE  6,45 

GO  10  996 

103  00  104  1*2,100 

IF  (TVI!)*  0.5*  M  -  X)  104,104,105 

104  CONTINUE 

107  WRirt  OUTPUT  TAPE  6,108 

108  FORMAT  (28H1  X  IS  GREATER  THAN  TV( 100 ) /1H0/IH0 ) 

GO  TO  996 

105  V*  VV(I-l) 

114  IF  (SENSE  SWITCH  1)  106,  110 

106  V*  V*  VO 

110  lFlX-TBIlM  111,112,112 

111  WRITE  OUTPUT  TAPE  6,55,X,TB(ll 
WRITE  OUTPUT  TAPE  6,45 

GO  TO  996 


t  * 

I 


f 


r 


C 

D 

0 


IF  (  X  -  TB(I)  -0.5*H  )  116,  111,  113 
113  CONTINUE 

WRITE  OUTPUT  TAPE  6,125 

125  FORMAT  I28H1  X  IS  GREATER  THAN  TB ( 100 1 / IH0/1H0 1 
GO  TO  996 

116  B  =  BB  (1-1) 

117  IF  (SENSE  SWITCH  41  118,  119 

118  B  *  B+BO 

119  WG  *  WG  -B*H 

IF  (WG)  109,  109,  120 
109  WRITE  OUTPUT  TAPE  6,15,WG 
PRINT  15, WG 
GO  TO  996 

120  IF  (  X  -  TCSI(l)  )  121,  l 22,  122 

121  WRITE  OUTPUT  TAPE  6,65, X.TCSI ( 1 1 
WRITE  OUTPUT  TAPE  6,45 

GO  TO  996 

122  DO  123  I*2,N$I 

I F ( X-TCS I  ( I  )  +  .5»H 1124,123, 123 

123  CONTINUE 

124  CSl*DSIDT( 1-1) •( X-TCS I( 1-1) )*CCSI< I-l) 

IF(X-TIR(  11  )126,127,127 

126  WRITE  OUTPUT  TAPE  6,81,X,T!RU) 

WRITE  OUTPUT  TAPE  6,45 

GOTO  996 

127  DO  128  I *2 ,NC !R 

IF(X-TIR(  I  )+.‘.*H  1129,128,  128 

128  CONTINUE 

129  CCIR-CIRC  t-i)*DCIRDT( I- 1 ) * ( X-TIR ( I- 1 ) ) 

132  I F l IY)130,153,130 

153  I Y* l 

GO  TO  1 40 

130  E  *  0. 

XS  «  X 

DO  131  1-1,6 

131  Y ( 1*10)  *  Y ( I ) 

CALL  RNGKTA  (H,  6,  1,  1) 

134  IF ( Y( 6) -V8H 1160,160,133 
160  PHI  -  V (6 * *•( 1./3. 1 

-  INCREMENTING  TOTAL  ENERGY  — — E2 - 

CALL  YPRiWE 

DK6«(WG+Y(51».5»RH0»Y(6n»Y(2l»F(2)2G 
DPE-  (RH0*Y(6)-WG-Y(5) )»Y<2) 

ODR AG-C0«.60375*PHI »PHI *KHO*Y IC)-A8SF(YI2)1»Y(2)/G 
0EWA-P»F(6) 

0EVE«Y(4)»CP*F(5) 

EP«EP»(-DK£tOPE-OORAG-OEWA4nEVE*OEH2-D6H34D£H4*OEM54DEH7)»M 

ET0T-CV»V(5l*Y(4)4Cl*WF0«Y(3) 

ENET-FTOT-ETOT 1 

ECHK-IENET-EPl/ENET 

FL-Y( 31 *( XMA»V (41/1 XMB»TA)-l • )-WG 

VALD»VALD*H»V/((XMA»Y(4))/(XM6»TA)-l.l 

GO  TO  140 


138  GO  TO  996 

133  DO  135  LAM8DA*1,  20 

E  *  £  ♦  <  V3M  -  Y  (  6  )  )  /  H 
X  =  XS 

00  136  1=1,6 

136  V ( I )  =  Y  ( I  + 10  ) 

CALL  RNGKTA  (  H,  6,  1,  1) 

IF(Y(6)-VBM)160»160»139 

139  IF  (SENSE  SWITCH  5)  137,  135 

137  PRINT  35  ,  X,  (  YU),  1  =  1,6),  V8M,  E 
135  CONTINUE 

145  PRINT  25 
PAUSE 

I F ( SENSE  SWITCH  5)  133,  138 

140  |FIY(1)-ZZI1)I141,142,14? 

141  IX  =  1 

GO  TO  143 

142  IF  I  X  -XT)  144,  141,  141 
144  I X=2 

IF(X-XP)157,143,143 

143  XP  *  XP  ♦  DPR 
IF(LP-17)146,147,147 

147  IF  (SENSE  SWITCH  3)  148,  149 

148  PRINT  994,  (TITLE  (I),  1*1,12) 

149  WRITE  OUTPUT  TAPE  6,  994,  (TITLE(I),  1-1, 12) 

WRITE  OUTPUT  1APE  6,  976 

LPC  *  LPC*1 
LP  *  0 

146  VR  *  V  -  RHO*E 

Ah* INTF (  ( GMT S*X )/ 3600*  1 

AM* I NTF ( (GMTS*X-3600.*AH)/60. ) 

AS*GMTS*X-60.*AM-3600.*AH 
AM*ABSF( AM) 

AS* A8SF ( AS ) 

PERB* 1 00. • ( Wl- WG ) /WT 

PERV*100.*VAL0/WH0 

PERE* 100. • ( WHO-Y (5))/(W80)-PFRV 

VEL*60.*Y(2> 

PM*0. 478801 •? 

IF  (  SENSE  SWITCH  3)  150,  151 

150  PRINT  9l,AH,AM,AS,Y(i),VEL,V,PERVfB,PER»,P6RE 

151  WRITE  OUTPUT  TAPE  6,977 , AH, AM, AS, X , Yl 1 ) ,V£L ,PN, TA , V( 3) , YI4) , Y(6 I »W 
lG,Y(S),fL.PERB,P€RV,PERfc,LAMB0A,6CHX 

LP*  LP* l 

155  I F f SENSE  SWITCH  2)156,157 

156  WRITE  OUTPUT  TAPE  6,5001 ,0EH2,0EH3,0EH4,DEH5,0EH6,0EH7 

157  CONTINUE 

IF(X-XPL)  210,209,209 
209  XPL*XPL*DP 
NPL*NPL*1 

IF(NPL-400)  211,211,210 
211  N»NPL 

Yl(N)*Y(l) 

V  31 N) *Y I  3 ) 

V4IN ) *V( 4 ) 


XXX ( N  )  *X 
TAA (N)*TA 

210  GOTO  {212,152), IX 
212  CALL  PLOT 
GOTO  906 

152  IF  (SENSE  SWITCH  6)  165,  101 

165  PRINT  166 

166  FORMAT  (22H0  RESET  SENSE  SWITCH  6/1H0/1HG) 
PAUSE 

996  READ  INPUT  TAPE  5, 10,  IT 
IF ( IT-2)  991,200,991 

991  1 F ( IT)  992,1001,992 

992  WRl’E  OUTPUT  TAPE  6,993 
CALL  EXIT 

300  READ  INPUT  TAPE  5,20, 1 T , I  XX ( I) , I* 1 , 7  ) 

I  P=3 

I F ( IT)  997,302,970 
302  Cl  =  XX( 1) 

C4=XX ( 2 ) 

C5=XX<3) 

C6=XX(4) 

ABIR=XX ( 5 ) 

EMI R=XX ( 6 ) 

ABUV=XX( 7) 

304  READ  INPUT  TAPE  >,20,  !T 
I F <  IT)  997,971,1000 

400  READ  INPUT  TAPE  5,  20,  IT,  I  XXII),  1=1,7  ) 
IP=4 

IF (  IT)  997,402,970 
402  WGO-XXU) 

WFO-XX  (2) 

VBM=XX( 3 ) 

FL0-XXI4) 

C0«XXC5! 

CB*XX (6;  * 

404  READ  INPUT  TAPE  5,  20,  IT 
.FUT)  997,971,1009 

500  READ  INPUT  TAPE  5,20, IT ,( XX (I) , 1*1 ,7 ) 

IP=5 

IF ( 3T  I  997,502,970 
502  7B0*XX(D 
TF0*XX(2) 

20*XX ( 3 ) 

Y20=XX(4) 

DEC* XX ( 5  ) 

Xf,AT*XX  (6) 

XLnNG=XX( 7 ) 

504  READ  IrtPUT  TAPE  5,  20,  IT 
l F (  IT)  997,971,1000 

950  READ  INPUT  TAPE  5,20, 1 T , I  XXII) , 1*1 ,7) 

I  P=6 

IF (  IT)  997,952,970 
95?  H  «  XX(1) 

DPR=XX ( 2 ) 

X0=XX 1 3 ) 


XT*XXI4) 
GHTS*XX 1 5) 
GHA=XX( 6) 
DP*XX 17) 


954 

REAO  INPUT  TAPE  5, 

I F  ( IT)  997,971,1000 

975, 

IT 

700 

IP=7 

DO  701  1*1,100 

READ  INPUT  TAPE  5, 
IF{  17)997,701,704 

30, 

IT, 

TCSI(I) 

,  ccsim 

701 

CONTINUE 

REAO  INPUT  TAPE  5, 
IF(IT)  997,971,706 

30, 

IT 

706 

1*1-  l 

704 

MS i = I - l 

DO  705  J=2 , Nb I 

US  IOT ( J-l )*( CCS! ( J ) 

-CCSKJ- 

■in/(TC$ItJ)-TCSl  (J 

705 

CONTINUE 

GOTO  1000 

800 

READ  INPUT  TAPE  5, 

I  P=8 

IF( 11)997,802,970 

30, 

IT, 

XX(  1) 

802 

VO  =  XXtl) 

DO  801  I*  1,100 

READ  INPUT  TAPE  5, 

I F  1  IT)  997,801,804 

30, 

IT, 

rvm. 

VVI I ) 

801 

CONTINUE 

READ  INPUT  TAPE  5, 

1 F (  IT)  997,971,805 

30, 

IT 

805 

1*1  +  1 

804 

NVI=I-l 

GO  TO  1000 

900 

READ  INPUT  TAPE  5, 

30, 

IT, 

XX(1  ) 

I  P=9 

I F (  IT!  997,902,970 

% 

902 

80  *  XX( 1) 

DO  901  1*1,100 

READ  INPUT  TAPE  5, 
IF (  I  T  )  997,901,905 

30, 

IT, 

TB( I ), 

BB  ( I ) 

901 

CONTINUE 

READ  INPUT  TAPE  5, 

I F  < IT)  997,971,904 

30, 

IT 

904 

1*1  +  1 

905 

NBI*I~1 

GO  TO  1000 

600 

READ  INPUT  TAPE  5, 
IP*10 

IHIT)  997,602,970 

30, 

IT, 

XX  (  1  ) 

602 

RHOO*  XX(l) 

DO  603  1*1,100 

REAO  INPUT  TAPE  5, 

30, 

IT, 

zzm, 

TZII) 

IF  (IT)  997,  603,  608 

603 

CONTINUE 

READ  INPUT  TAPE  5, 

30, 

IT 

IF  I  I T  J  997,971,609 
609  I  *  I  +  l 

608  Nr  *  i~i 

RTZUI  *  RH00  •  TZ €  1 » 

00  604  1*2 ,  NT 

IF  (  TZII)  -TZtI-lM  605,  606,  605 

605  TZA(I-i)  =  (  TZII)  -  TZU-m/IZZU  >  -ZZII-1H 

TZ9(i-n  *izzm  •tzii-1)  -  zzu-n  *  rzun  /  izzm  -  zzu-m 

RTZtn  *  RTZ(l-l)  «ITZU-l)/TZm  l**<  XMA/(  R  •TZAU-l))  ) 
GO  TO  604 

606  rZAI  l-l)  *  0.0 

T ZB (  I-n  -  TZIII 

RTZ  ( I  )  *  RTZU-1I  +EXPF  I  XMA«  ( ZZ  I I-U-ZZUI 1/ (R  *TZII-I)n 
604  CONTINUE 
GO  TO  1000 

960  IP*U 

DO  961  1*1,100 

READ  INPUT  TAPE  5, 30,  I  T  »Z  IR  ( I)  ,RT  IR  tii 
FIR  {I  i*SBOLZ*RT IR( I )*»4.0 
IF! IT) 99 7, 96 1,962 

961  CONTINUE 

READ  INPUT  TAPE  5,30, IT 
I F I  IT)  997,971,964 
964  1*1+1 

962  NIR® l~l 

DO  963  J*2,NIR 

DIRDZ  { J-l  )*  (FIR1  J)-FIRIJ-in/IZIR(  JI-ZIRIJ-in 

963  CONTINUE 
GOTO  1000 

980  IP* 12 

DO  984  MZY*i„400 
ET!MF(MZY)*0. 

EALT(MZYI*0. 

984  CONTINUE 

DO  981  1*1,400 

REAO  INPUT  TAPE  5,31 , 1  T.ETIMEI  I)  ,EAITU) 

I F I  IT) 99 7, 981, 982 

981  CONTINUE 

READ  INPUT  TAPE  5,30,  IT 
IF  I  IT)  997,971,983 
983  1*1+1 

982  NEX*I-l 

GO  TO  1000 

1400  I P* 1 3 

DO  1401  1*1,100 

REAO  INPUT  TAPE  5,30,  I T  ,  T  IR  (  I)  ,  CIRU) 

I F ( ITJ997, 1401, 1402 

1401  CONTINUE 

REAO  INPUT  TAPE  5,30, IT 
IFUTI  997,971,1403 

1403  1*1+1 

1402  NC IR* I~l 
001404  1*?,NCIR 

0C!RDTU-l)*fCIR(I  )-CIR(I-m/{TIRI  II-TIRI  1-1)1 

1404  CONTINUE 
GO  TO  1000 


SUBROUTINE  RNGXTA { HI ,N1 , N2 »N3 ) 
COMMON  X,Y«  F,Q 

DIMENSION  YlIOOltFIIOO),  G(IQO) 
IFiN3-l)2»l,2 

1  H*H1 
HH=.5*H 
N*N1 
M=N2 

DO  3  1=1, N 
3  8(11*0,0 

2  DO  4  J*1,M 
CALL  YPRIME 
DO  5  1=1, N 
S=F( I  )«H 

T  =  .5*IS-2.*Om  I 
Y(  1 1 =Y 1 1 ) +T 

5  6!I)*0I I »+3.«T-.5*S 
X=X+HH 

CALL  YPRIME 
DO  6  1*1,  N 
S=F  ( I  )*H 

T=«2 9289322*1 S-QI I )  ) 

Y { I )*YI I )+T 

6  QII)*0(l)O.iT  -.29289322*S 
CALL  YPRIME 

DO  7  1*1, N 
S=F(n*H 

T*1.7071067*(S~0m  ) 

Y ( I l=Y( I»*T 

7  Q(  I)  =  0m  +  3.*T~1.707106*S 
X=X+HH 

CALL  YPRIME 
DO  8  1*1, N 
S=F ( I ) *H 
T=(S-2.*Q( I ))/6. 

Y ( I ) *Y( I  )  +T 

8  Q( I )=Q( II+3.*T-,5*S 
A  CONTINUE 

RETURN 

END 


SUBROUTINE  YPRIME 
COMPILATION  80 

10  FORMAT ( 14H  DIVIDE  CHECK  ,111 
DIMENSION  SCALE (60) 

DIMENSION  ZZ(IOO),  TZ(IOO),  RTZUGO),  TZA(IOO),  TZB(IOO) 
DIMENSION  TCSI(100),CCSI(100),  TVUOO),  VV(IOO) 

DIMENSION  TB(IOO),  8B(100)  ,  XX(IQO) 

DIMENSION  Y(100),  F(IOO),  C(IOO) 

DIMENSION  TITLE  (12) 

DIMENSION  ZIR(100),FIR(100),DIRDZ(100),DSIDT(100) 

DIMENSION  ZUV(  100) ,FUV( 100) ,DUVDZ ( 100) 

DIMENSION  Yl( 400) , Y3 1 400 ) , Y4( 400 > , XXX( 400 ) » TAA (400) 

DIMENSION  £  T IME  ( 400  ) »  EAL  T  (  400 ) 

COMMON  X,  Y,  F,  0 

COMMON  NT,  11,  TZ,  RTZ,  TZA,  TZB 

COMMON  Cl,  C2,  C3,  C4,  C5,  C6,  C7 

COMMON  TA,  RT,  P,  RHO,  PHI,  WG,  OMEGA,  CSI,  E,  V 

COMMON  XMA,  XMB ,  G,  CP,  CV,  R,  CD,  VBM  ,WF 

COMMON  DEH2 ,DEH3, DEH4, 0EH5 » DEH6 » DEH  7 

COMMON  WFH,DWFH, WFO 

COMMON  CB 

COMMON  ZIR,FIR,DIRDZ,NIR,ZUV,FUV,DUVDZ,NUV,DSIDT 
COMMON  ABIR,EMIR, ABUV 
COMMON  M,Y1,Y3,Y4,TAA,XXX 
COMMON  DP,XT,XO,ETIME,EALT 
COMMON  CCIR 

COMMON  SLATO,CLATD,RAD,GHA,XLONG, PZER, FSOL, SBOLZ, AH,AM,AS,GMTS 
CALL  RHOT  (  Y(l),  TA,  RT  ) 

P=RT*K/ XMA 

RHO  =  RT  /  TA 

PHI  *Y ( 6 ) •* ( 1,0/3. 0) 

DO  181  1=2, NIR 

IF(Y(1)-ZIR(I  )) 182, 181, 181 

181  CONTINUE 

GO  TO  ( 184,183),IXY 
184  PRINT  20, X, Y( 1 ) 

20  FORMAT ( IHO, 9H  AT  T I ME*F7. 0, 75H  SEC  THE  BALLOON  KENT  ABOVE  THE  INFR 
IARED  FIELD  SPECIFIED  IN  THE  !NPUT  DATA./18H  BALLOON  ALTI TUDE-F7.0, 
26H  FEET./123H  START  TO  CONTINUE  FLIGHT  USING  INFRARED  DATA  OF  UPPE 
3R  POINT  UF  SPECIFIED  FIELD  FOR  AS  LONG  AS  BALLOON  REMAINS  ABOVE  FI 
4ELD./51H  SENC6  SWITCH  6  DOWN  AND  START  TO  TERMINATE  FLIGHT.) 

PAUSE 
IXY«2 
GO  TO  163 

182  IXY-l 

183  FLUX  I R»  DIRDZ( !-l)*!Y( 1 )-Z!R( I-l> )*FIR( l-ll 
DEH2»CCIR*ABIR*4.B3»PHI#PHI*FLUXIR 
DEH3*EHIR»4. 83«PHl »PHI *SBQLZ  ® ( Y I 3)**4. I 

186  XLAH«RAD* (GHA-XL0NG+X/240. ) 

CAM*  SLATD  ♦  CLATO*  COSF(XLAH) 

IF  «CAM)  100,101,101 

101  AIRM*  CP/PZER)*II 1228.6  4376750. 44*CAM*CAM)«*. 5  -613.8«CAM) 

TRANS*. 5*(EXPF(-.65»AIRM)4EXPF(-.095*AIRM)) 

FLUXUV-  FSOL»TRANS 


GOTO  no 

100  RETH  *  20903*520* 

CAL*RETH/ ( RE  TH+Y (  1 )  ) 

CAMM—l.M i.~CAL*CAL )  **.  5 
IF  !CAMM-CAM)102» 103,103 

102  ZZ 1  = I RETH+Y I 1 ) ) *  1 1 .-CAM*C AM ) *• . 5-RE TH 
CALL  RHOT (ZZi»TTT»RTT) 

PZl=RTT*R/XMA 

AIRMl=35.1*PZl/PZER 

AIRM2= I P/PZER ) • I I 1228.6+376750.44*CAM*CAM)»».5*613.8*CAM) 
TRANI=.5*(EXPF{-.65«AIRM1)+EXPF(-.095«AIRM1 ) ) 

TRAN?=.5*IEXPF l-.65*AIRM2)+EXPF(-.095»AIRM2 ) ) 

FLUXUV=FSOL»  TRAN1 *TR AN 1/TRAN2 
GOTO  110 

103  FLUXUV-O. 

110  CONTINUE 
0EHH7=DEH7 

DEH7=CSI*ABUV*1.208*PHl#PHI»FLUXUV 

IF(DEHH7)iU,ill,il5 

111  IF(DEH7)il9»ll9,ll4 

114  IF(X)U9,  119,120 

120  INM*1 

GO  TO  121 

112  WRITE  OU I  PUT  TAPE  6, 11 3 , AH, AM, AS 

113  FORMAT  I 1H0»?5X,4HGMT *3F3*0,5X» 18HSUNRI SE  AT  BALLOON) 

GO  TO  119 

115  1F(DEH7)118, 118,119 

118  I NM=2 

121  AH* I NTF I ( GMTS+X ) /3600. ) 

AM-INTFI (GMTS+X-3600.*AH>/60. ) 

AS-GMTS+X-60. • AM-3600. • AH 
AM-ABSFI AM) 

AS-ABSF ( AS ) 

GO  T0I112, 116), INM 

116  WRITE  OUTPUT  TAPE  6, 117, AH, AM, AS 

117  FORMAT  I 1H0,25X,4HGMT-3F3.0»5X» l 7H SUNSET  AT  BALLOON) 

119  VISA*l.096E-05*(TA/460. )**.883 
CA|R»260.l«VISA 

VlSH*l.2llE-05*<Y(4)/460. )»».682 
CHE-1443. <*VISH 
RHQH«<XMB*P)/(R»Y<4) ) 

FNUS».37»( I l.24*RH0«PHl«ABSFIYl?) )/VISA)«».6) 
DEH4«3.895*C4»CAIR»PHl*FNUSMTA-Yf  3>) 

XNUSA-. 1612*PHI *( ( RHO»RHO»G*.7»ABSF I TA-Yt  3 ) ) )/l VISA»VISA*TA) >** 
l( l .0/3.0) 

0EH5- 3.895»C5*CA!R»PHI*XNUSA*ITA-Y(3)) 

XNU$H«.16i2*PHl*URHOH»RHOH*G*.7«ABSFtY(4)-Y<3) )  )/(VlSH»VISH*Y(4) ) 
l)«*(l. 0/3.0) 

0EH6*3.895*C6*CHE»PHI*XNUSH»tY<4)-YI3l) 

Fit)  *  Y(2) 

F 12 ) - I IRH0*YI6) -WG~Y| 5) l•G-.6038*C0•PH^•PHI *RHO*Y( 2 )*AB$F IY<2) ) )/ 

1 I WG+Y ( 5 )4CB*RH0»V 1611 

FC3)»I 06H2-OEH340EH440EH540EH640EH7 )/l Z UWF ) 

F 1 4 )- ( -DEH6-RH0*YC 6)*YI2))/(V(5 )* I CV+R/XMfi) ) 

F 1 5 1 * RH0H»E-AH0H* V/ 1 RHO-RHOH ) 


F(6I*IIFC4)/Y(4))4|F(4)/Y|5) )♦ ( XWA*Y (2 ) /( R»  TA) > |*Y(6> 

RETURN 

END 


SUBROUTINE  RHOT  (  Z,  T,  RT) 

C  COMPILATION  80 

DIMENSION  SCALE ( 60 ) 

DIMENSION  Z2  I  100)  ,  TZ(IOO)*  RTZ(IOO),  TZA(IOO),  TZB(IOO) 
DIMENSION  TCSI(100),CCSH100),  TV(IOO),  VV(IOO) 

DIMENSION  TB(IOO),  88(100)  ,  XX(IOO) 

DIMENSION  Y ( 100 } »  F(IOO),  0(100) 

DIMENSION  TITLE  (l?) 

DIMENSION  ZIR( 100)tFIR( 100) tDIROZ(lOO) tOSIDT(lOO) 

DIMENSION  ZUV ( 1 00 ) , FUV ( 100 ) , DUVDZ (100) 

DIMENSION  Y1 ( 400 ) » Y3( 400) , Y4( 400 ) , XXX ( 400 ) *  TAA ( 400 1 
DIMENSION  ETIME(400),EALT(400) 

COMMON  X,  Y,  F,  Q 

COMMON  NT,  lit  TZ»  RTZ,  TZA,  TZB 

COMMON  Cl,  C2,  C3,  C4,  C5,  C6,  C7 

COMMON  TA,  RT,  P,  RHO,  PHI,  WG,  OMEGA,  CSI,  E,  V 

COMMON  XMA,  XMB ,  G,  CP,  CV,  R,  CD,  VBM  ,WF 

COMMON  DEH2 ,DEH3,DEH4, DEH5,DEH6,DEH7 

COMMON  WFH, OWFH, WFO 

COMMON  CB 

COMMON  21R, FIR, DIRD2,NI R, ZUV, FUV, DUVDZ «NUV,DS IDT 
COMMON  ABIR,EMIR,ABUV 
COMMON  M,Yl,Y3,Y4,TAA,XXX 
COMMON  DP,XT,XO,ETIME,EALr 
IF  (  Z-  ZZ ( 1 )  )  100-200,  200 
100  T  *  TZ ( l ) 

RT  *  RTZ(  1) 

RETURN 

200  NMl  *  NT-l 

DO  201  J*l,  NMl 

IF  (  Z  -  ZZC  J*l»)  300,  300,  201 

201  CONTINUE 

GO  TO  ( 202,203 ) , I YX 

202  PRINT  20, X, Y( 1 ) 

20  FORMAT ( IH0,9H  AT  TIME*F7.0,63H  SEC  THF  BALLOON  WENT  ABOVE  ATMOSPHE 
IRE  SPECIFIED  IN  INPUT  DATA. /IBM  BALLOON  ALT  I TUDE»F7.0,6H  FEET./128 
2H  START  TO  CONTINUE  FLIGHT  USING  ATMOSPHERIC  DATA  OF  UPPER  POINT  0 
3F  SPECIFIED  ATMOSPHERE  FOR  AS  LONG  AS  BALLOON  REMAINS  ABOVE  IT./51 
4H  SENSE  SWITCH  6  COWN  AND  START  TO  TERMINATE  FLIGHT.) 

I  YX*2 
PAUSE 
GO  TO  203 

300  I YX« l 

203  IF  (  TZA(J))  301,  302,  301 

301  T  -  TZ A(  J)  *  Z  ♦  TZB  IJ) 

RT*RTZ( J)* I TZ ( J)/ Tl XMA/(R*TZA( J ) ) ) 

RETURN 

302  T  «  TZB  (J) 

RT«RTZ( J)»EXPF ( XMA*C ZZ ( J )-Z )/( R*T ) ) 

RETURN 

ENO 


SUBROUTINE  PLOT 


SUBROUTINE  PLUT 
C  COMPILATION  70 
3000  FOR MAT ( 14H  GRAPH  PRINTED) 

DIMENSION  SCALE  16'") 

DIMENSION  iznc‘'),n(  1 00 ) • R  T2 ( 1 00 ) *  T  Z  A ( 1  OP ) *  TZ  0 ! I DC ) 
DIMENSION  TCSH  lnP),CCSI(  1 00  5  *  TV  I  100)  fVV(iO0) 
DIMENSION  TCI  10'  )  ,BH(  100)  ,  XX  ( 1  Of') 

DIMENSION  Y(  ICO)  «F  I)  00.)  ,0(100) 

DIMENSION  TITLE  (12) 

DIMENSION  ZIRHOO)  ,riR(  IOO),DIKOZ(  100)  ,  OS  IDT  (100) 

D I  MENS 1  ON  ZUV { 1 00 ) , Ft  I V ( 100 ) , nuvnz ( 1 : ' ) 

DIMENSION  Y1  (400)  ,Y3(400)  ,Y4(40t  »  ,XXX(4<:C) »  TAA  (400 
0 1  MENS  I UN  ET I  ME ( 400 ) , t  ALT ( 40D ) 

COMMON  *,Y,F,P 

COMMON  NT.ZZf TZ,RTZ,TZA,TZR 

COMMON  Cl  ,C?,C?»C4»C'.>»Ct,C7 

COMMON  TA,RT, P, RHU, PH  I , WG, OMEGA , C S  I  ,c, V 

COMMON  XMA,  XMK,  G, CP,  C  V,  R ,  CD,  Vt'.M ,  WF 

COMMON  CEH2,0L‘H3,PgH4,DFH5,DEM6,DEH7 

COMMON  VlFH,  f>WFH»WFO 

COMMON  CD 

COMMON  ZIR,FIR,niRO?,NlR,  ZUV, FUV, I'UVOZ , NUV»OSI PT 

COMMON  AB IR v £M IP , ARUV 

COMMON  M,Y1,Y3,Y4,TAA,XXX 

COMMON  t)P,XT,XO,CTIME,EALT 

SC  ALE  I  5 ) =6H  325 

SCAL  E  (  7 )  =  4H  more 

SCALE(ll)=6H  35C 

SCALE (13 )  =6M  2000"' 

SC  ALE ( l 7 ) =6H  375 
SCALE ( 19 ) «6H  30000 
SCALE (23) *6H  40C 
SCALE ( 25 ) *6H  4000^ 

SCALE!2P)*6H  425 
SCALE (ii)*6H  50000 
SCALE. i 37 )*6M  450 
SCALE! 37J *6H  60000 
SCALE ( 41 »*6H  475 
SC  ALE (4  31 *6H  7CC0C 
SCALEt47)-*6H  50C 
SC  ALE  1 49 1 »6H  9C000 
SCALE ! 53) «6M  5?5 
SCALEC5>)«6H  90000 
SCALE! 59)*6H  550 
16  M»1 

MT»XINTF I (XT-XO)/IOHOQ* )♦! 

7  MM«M*1 


GO  TO  (4, l? , 13, 14, 15, 16,1?) *HM 
12  SCALE!  1  )*6H  0-3 


GO  10  9 

13  SCALE ! 1  )*6H  3-6 
GO  TO  9 

14  SCALE ( 1  )«6H  6-9 
GO  TO  9 

15  SCALE! 1 )*4H  9-12 
GO  TO  9 


SUBROUTINE  PLOT 


16  SCALE(i‘=6H  12-15 
GO  TO  9 

IT  SC  ALE  1 1 ) =6H  15-18 
9  XM=FLOATF I M ) 

XXM=I XM-1 . )*10800. 

XXXM*XM*1 0800. 

CALL  LIMITS  IXXM,XXXM,0., 100000.) 
CALL  GRID  (XXM, 1800., 0., 10000. ) 

K=  XF l XF  I  (  1.0800./0P)*t  XM-1.  )  )♦! 

L  =  XF1XFI  ( 10800. /DP )*( XM) )+l 
DC30J=K,L 

90  CALL  POINTS  ( XXXC J) ,Y1 ( J) , 35 )  , 

1C  CALL  LIMITS  ( XXM, XXXM,  500. ,550.) 
00100J=K,L 

100  CALL  POINTS  I  XXXI J) ,Y3( J) , 15) 

11  CALL  LIMITS  ( XXM, X XXM, 300 *,550.  ) 

DO  110 J=K , L 

110  CALL  POINTS  I XXXI J ) ,Y4( J) , 1 1 ) 

CALL  LIMITS  l XXM, X XXM, 300 • , 550 •  ) 

DO 1?0 J=K, L 

120  "ALL  POINTS  IXXXI J),TAAI J),10) 

CALI  LIMITS  IXXM,XXXM,0., 100000.) 
OC130J* 1 , 400 

130  CALL  POINTS  ( ET IME I J) . EAL T I J ) , 14) 
CALL  GR  I  SCALE) 

4  PRINT  300 '} 

M*M+1 

IF  I M- 7)131, 19, ’9 

131  IF  IM-NT ) 7, 7, 19 
19  RETURN 

tNCIl»0»0»0»l, 0*0,0, 0,0,0»C,0, 0,0) 


SUBROUTINE  PLOT 

C  COMPILATION  FOR  147500  FT 
3000  FORMAT! I4H  GRAPH  PRINTED) 

DIMENSION  SCALE ( 60 ) 

DIMENSION  22  ( 100 )  t  TZ(IOO),  RTZ(IOO),  TZA(IOO)*  TZBUOO) 
DIMENSION  TCSI!100),CCSU100),  TV(IOO),  VV(IOO) 
DIMENSION  TB(IOO),  BB( 100)  »  XXtlOO) 

DIMENSION  Y( 100 ) »  F(100),  Q(IOO) 

DIMENSION  TITLE  (12) 

DIMENSION  2IR11 00 )»FIR(1 00 )fDIR 02(100) ,DSIDT( 100) 

DIMENSION  2UV! 100 ) * FUV( 100 ) *DUVD2 ! 100) 

DIMENSION  Y l ( 400 ) » Y3( 400) » Y4(400 ) *  XXX! 400 ) « TAA (400) 
DIMENSION  6T IME ( 400 ) « EALT ( 400) 

COMMON  X#  Y,  F,  0 

COMMON  NT,  2Z ,  TZ,  RTZ,  TZA,  T2B 

COMMON  Cl,  C2 ,  C3,  C4,  C5,  C6,  C7 

COMMON  TA,  KT,  P,  RHO,  PHI,  HG,  OMEGA,  CSI,  E,  V 

COMMON  XMA,  XM8 ,  G,  CP,  CV,  R,  CD,  VBM  ,WF 

COMMON  DEH2,DEH3,0EH4, DEH5«DEH6«DEH7 

COMMON  WFH, OWFH, WFO 

COMMON  CB 

COMMON  2 1R , FIR ,DIR02 ,N IR, 2UV,FUV,DUVD2 , NUV,DS!DT 
COMMON  ABIR,EMIR, A8UV 
COMMON  M,Y1,Y3,Y4,TAA,XXX 
COMMON  DP,XT,XO,ETIME,EALT 
SCALE ( 8 )*6H  340 
SCALE ( 9  )*6H  20000 
SCALE ( 16)»6H  380 
SCALE! 17 )»6H  40000 
SCALE! 24 )*6H  420 
SCALE !25)«6H  60000 
SCALE ( 32 )*6H  460 
SCALE! 33) -6H  80000 
SCALE!40)*6H  500 
SCALE (41 )*6H 100000 
SCALE! 48  >*6H  540 
SCALF(49)*6H120000 
SCALEC56»*6H  580 
SCALE ( 57 )*6H140000 
18  M*1 

MTaXlNTF ( (XT-XO)/ 10800* )♦! 

7  MM»M+1 

GO  TO  (4, 12,13, 14,15, 16, 17), MM 

12  SCALE! 1 )*6H  0-3 
GO  TO  9 

13  SCALE ! 1 )*6H  3-6 
GO  TO  9 

14  SCALE! 1 )*6M  6-9 
GO  TO  9 

15  SCALE 1 1 )«6H  9-12 
GO  TO  9 

16  SCALE! 1 1»6H  12-15 
GO  TO  9 

17  SCALE! 1 )*6H  15-18 


9  XM«FLOATF(M) 

XXM«  ( XM-1 . ) • 1 0800. 

XXXM»XM*10800. 

CALL  LIMITS(XXM,XXXM*0.  *147500* ) 
CALL  GRID  ( XXM, 1800., 0., 20000. ) 
K*XFIXF( (10800. /DP)* ( XM-i. ) )♦! 
L*XF I XF ( ( 10800. /DP )*( X# ) )*1 
DO  90  J»K,L 

90  CALL  POINTS  ( XXX( J ) ,Y l ( j ) , 35) 

10  CALL  LIMITS  ( XXM.XXXM, 300. * 595. ) 
DO  ICO  J*K  ,L 

100  CALL  POINTS  ( XXX( J ) ,Y3( J) , 15) 

11  CALL  LIMITS  ( XXM*XXXM, 300. ,595. ) 
011  110  J«K,L 

110  CALL  POINTS  ( XXX( J ) * Y4 ( J )  *  1 1 ) 

CALL  LIMITS  ( XXM,XXXM, 300. , 595. ) 
DO  120  J»K,L 

120  CALL  POINTS  ( XXX( J ) *TAA( J ) , 10) 
CALL  LIMITS  (XXM, XXXM,0., 147500. ) 
00  1310  J*1 ,400 

130  CALL  POINTS  ( FT IME (J ) , EALT ( J) , 14 ) 
CALL  GRAPH  (SCALE) 

4  PRINT  3000 
M*M*1 

IFIM-7) 131,19,19 

131  IF ( M-MT ) 7, 7, 19 
19  RETURN 

ENO 


